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Abstract: In these lectures I present a review of non-perturbative instanton effects in quantum 
theories, with a focus on large N gauge theories and matrix models. I first consider the structure 
of these effects in the case of ordinary differential equations, which provide a model for more 
complicated theories, and I introduce in a pedagogical way some technology from resurgent 
analysis, like trans-series and the resurgent version of the Stokes phenomenon. After reviewing 
instanton effects in quantum mechanics and quantum field theory, I address general aspects 
of large N instantons and then present a detailed review of non-perturbative effects in matrix 
models. Finally, I also consider two applications of these techniques in string theory. 
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1. Introduction 



Many series appearing in Physics and in Mathematics are not convergent. For example, most 
of the series obtained by perturbative methods in Quantum Field Theory (QFT) turn out to be 
asymptotic, rather than convergent. The asymptotic character of these series is typically an indi- 
cation that non-perturbative effects should be "added" in some way to the original perturbative 
series. 

The purpose of these lectures is to provide an introduction to asymptotic series and non- 
perturbative effects. Since this is already quite a wide topic, we will restrict our discussion in 
various ways. First or all, we will consider non-perturbative effects of the instanton type, i.e. 
effects due to extra saddle points in the path integrals (in particular, we will not consider effects 
of the renormalon type). Secondly, we will be particularly interested in the interaction between 
non-perturbative effects and large N expansions. Finally, in our discussion we will rely heavily 
on "toy models" of non-perturbative effects, and in particular on matrix models. This is not as 
restrictive as one could think, since matrix models underlie many interesting quantities in string 
theory and supersymmetric QFTs. 

A typical asymptotic series in a coupling constant g s , with instanton corrections, has the 
following form, 

£ a n g n s + e~ A ^ £ a«<£ + 0(e" 2 ^)- (1-1) 

n n 

Here, the first sum is the original divergent series. The second term is an one-instanton contri- 
bution. It has an overall, non-pertubative exponential in g s , multiplying another series (which 
in general is also asymptotic). The third term indicates higher instanton contributions. In order 
to understand this type of quantities, we will develop a three-step approach: 

1. Formal: we want to be able to compute the terms in the original, perturbative series, as 
well as the "non-perturbative quantities" characterizing the instanton corrections. These 
include the instanton action A, as well as the series multiplying the exponentially small 
terms. Notice that the resulting object is a series with two small parameters, g s and e^ A ^ 9s , 
which should be regarded as independent. Such series are called trans-series. Therefore, 
the first step in the understanding of non-perturbative effects is the formal calculation of 
trans-series. 

2. Classical asymptotics: The series (1.1) above is a purely formal expression, since already the 
first series (the perturbative one) is divergent. One way of making sense of this perturbative 
series is by regarding it as an asymptotic expansion of a well-defined function. Finding such 
a representation, once the original function is given, is the problem which I will call "classical 
asymptotics." In classical asymptotics, exponential corrections are ill-defined and are never 
written down explicitly, but they are lurking in the background: indeed, it might happen 
that, as we move in the complex plane of the coupling constant, an instanton correction 
which used to be exponentially small becomes of order 1. This is the basis of the Stokes 
phenomenon. 

3. Beyond classical asymptotics: once the classical asymptotics is understood, one can try to go 
beyond it and use the full information in the formal trans-series to reconstruct the original 
function exactly. A general technique to do this is Borel resummation. The combination of 
the general theory of trans-series with Borel resummation gives the theory of resurgence of 
Jean Ecalle, which is the most powerful method to address this whole circle of questions. 
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In these lectures we will follow this three-step program in problems with increasing order of 
complexity, from simple ones, like ordinary differential equations (ODEs, to more difficult ones, 
like matrix models and string theory. Let us quickly review these different problems in relation 
to the approach sketched above: 

1. In the theory of ODEs with irregular singular points, formal solutions are typically divergent 
series. In this case, the three steps are well understood. The first step is almost automatic, 
since the terms in the trans-series can be computed recursively, including the instanton 
corrections. The understanding of the classical asymptotics of these solutions is a venerable 
subject, going back to Stokes. The treatment beyond classical asymptotics is more recent, 
but has been well established in the work of Ecalle, Kruskal, Costin, and others, and we will 
review some of this work here. A closely related example is the calculation of integrals by 
the method of steepest descent. Here, the different trans-series are given by the asymptotic 
expansion of the integrals around the different critical points. The classical asymptotics 
of these series is also a well studied subject, and the resurgent analysis has been also 
considered by Berry and others. 

2. In Quantum Mechanics (QM) and Quantum Field Theory, the formal procedure to generate 
the leading asymptotics series is simply standard perturbation theory. Instanton corrections 
are obtained by identifying saddle-points of the classical action, and by doing perturbation 
theory around this instanton background one obtains exponentially small corrections and 
their associated series. The study of classical asymptotics and beyond is more difficult, 
although many results have been obtained in QM. In realistic QFTs, perturbation theory 
is so wild that it is not feasible to pursue the program, but in some special QFTs -namely, 
those without renormalons, like Chern-Simons (CS) theory in 3d or N = 4 Yang-Mills 
theory- there are some partial results. In these lectures I will consider some aspects of the 
problem for CS theory. 

3. An interesting, increasing level of complexity appears when we consider large N gauge 
theories. Here, the computation of formal trans-series becomes more complicated, since 
in the 1/N expansion we typically have to resum an infinite number of diagrams at each 
order. For example, in standard QFT, the instanton action is obtained by finding a classical 
solution of the equations of motion (EOM) with finite action, while the action of a large 
N instanton is given by the sum of an infinite number of diagrams. Another way of 
understanding this new level of complexity is simply that in large N theories we have two 
parameters in the game, N and the coupling constant g s , or equivalently g s and the 't Hooft 
parameter t = Ng s . Due to this, the Stokes phenomenon of classical asymptotics becomes 
more complicated and leads to large N phase transitions (in fact, although it is not widely 
appreciated, standard phase transitions are examples of the Stokes phenomenon; see [107] 
for a nice discussion of this). However, in the simple toy case of large N matrix models, 
the 1/N expansion is still very close to the ordinary asymptotic evaluation of integrals, and 
we will be able to offer a rather detailed picture of non-perturbative effects. 

4. Finally, in string theory things are even more complicated. Even at the formal level we 
cannot go very far: there are rules to compute the genus expansion, but the rules to do 
(spacetime) instanton calculations are rather ad hoc. We know that in general instanton 
corrections involve D-branes or p-branes, and we know how to obtain some qualitative 
features of their behavior, but a precise framework is still missing due to the lack of a 
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non-perturbative defition. An alternative avenue is to use large N dualities, relating string 
theories to gauge theories, in order to deduce some of these effects from the large N duals. 
This makes possible to compute formal trans-series for non-critical strings and some simple 
string models. 



The plan of these lectures is the following: in section 2, I will review some aspects of asymp- 
totics series and in particular of the series appearing in the context of ODEs. This includes 
formal trans-series, Borel resummation, the Stokes phenomenon, and the connection between 
large order behavior and trans-series. I have tried to provide as well some very elementary ideas 
of the theory of resurgence. In section 3, I review non-perturbative effects in QM and in QFT. 
The results in QM are well known to the expert; they have been analyzed in detail and made 
rigorous in for example [46] . The study of non-perturbative effects in QFT from the point of view 
advocated in these lectures is still in its infancy, and I have contented myself with explaining 
some of the results for CS theory, which are probably not so well known. In section 3 I also 
introduce some general aspects of instanton effects in large N theories which are probably well 
known to experts, but difficult to find in the literature, and I illustrate them in the simple case 
of matrix quantum mechanics. Section 4 is devoted to the study of non-perturbative effects in 
matrix models, starting from the pioneering works of F. David in [44] and explaining as well 
more recent results where I have been involved. In section 5, I present two applications of the 
techniques of section 4 to string theory, by using large N dualities. Finally, in section 6 I make 
some concluding remarks. 



2. Asymptotics, non-perturbative effects, and differential equations 

2.1 Asymptotic series and exponentially small corrections 

A series of the form 

oo 

= (2.1) 



n=0 



is asymptotic to the function f(w), in the sense of Poincare, if, for every N, the remainder after 
N + 1 terms of the series is much smaller than the last retained term as w — > 0. More precisely, 



lim w N 

w—>0 



' N \ 

f(w)-J2^w n =0, (2.2) 
v, n=0 / 



for all N > 0. Notice that, in an asymptotic series, the remainder does not necessarily go to zero 
as N — > oo for a fixed w, in contrast to what happens in convergent series. Analytic functions 
might have asymptotic expansions. For example, the Stirling series for the Gamma function 

is an asymptotic series for \z\ — > oo. Notice that different functions may have the same asymptotic 
expansion, since 

f(w) + Ce~ A/w (2.4) 
has the same expansion around w = than f(w), for any C, A. 
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In practice, asymptotic expansions are characterized by the fact that, as we vary N, the 
partial sums 

N 

S N (w) = Y,a n w n (2.5) 

n=l 

will first approach the true value f{w), and then, for N sufficiently big, they will diverge. A 
natural question is then to obtain the partial sum which gives the best possible estimate of f(w). 
To do this, one has to find the N that truncates the asymptotic expansion in an optimal way. 
This procedure is called optimal truncation. Usually, the way to find the optimal value of N is 
to retain terms up to the smallest term in the series, discarding all terms of higher degree. Let 
us assume (as it is the case in all interesting examples) that the coefficients a n in (2.1) grow 
factorially at large n, 

a n ~ A~ n n\ . (2.6) 
The smallest term in the series, for a fixed \w\, is obtained by minimizing N in 



N 

\a n w n \ = cN\ 



(2.7) 



By using the Stirling approximation, we rewrite this as 

cexp|iV(log JV-l-logX)}, (2.8) 

where 

X 



A 

w 

The above function has a saddle at large N given by 



(2.9) 



N* = X. (2.10) 

If\w\ is small, the optimal truncation can be performed at large values of N, but as \w\ increases, 
less and less terms of the series can be used. We can now estimate the error made in the optimal 
truncation by evaluating the next term in the asymptotics, 

e( w ) = C N , +1 \w\ N ' +1 ~ e-l^l. (2.11) 

Therefore, the maximal "resolution" we can expect when we reconstruct a function f(w) from 
an asymptotic expansion is of order e(w). This ambiguity is sometimes called a non-perturbative 
ambiguity. The reason for this name is that perturbative series are often asymptotic, therefore 
they do not determine by themselves the function f(w), and some additional, non-perturbative 
information is required. Notice that the absolute value of A gives the "strenght" of this ambiguity. 

It is instructive to see optimal truncation at work in a simple example. Let us consider the 
quartic integral 

1(g) = dze" z2 / 2 -^ 4 / 4 , (2.12) 

which is well-defined as long as Re(g) > 0. The asymptotic expansion of this integral for small g 
can be obtained simply by first expanding the quartic term in the exponent and then integrating 
term by term. One obtains, 



oo 



I(g) = J2^9 k , (2.13) 

fc=0 
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Figure 1: We illustrate the method of optimal truncation for the quartic integral (2.12) by plotting the 
difference (2.16) between the integral and the partial sum of order N of its asymptotic expansion, as a 
function of N, for g = 0.02 (left) and g = 0.05 (right). 



where 

(_ 4 )-fc r°° z ik ^_ fc (4A;-l)!! 

This series has a zero radius of convergence and provides an asymptotic expansion of the integral 
1(g). The asymptotic behavior of the coefficients at large k is obtained immediately from Stirling's 
formula 

a k ~ (-A) k k\, (2.15) 
therefore \A\ = 1/4 in this case. In Fig. 1 we plot the difference 



\I(g)-S N (g)\ (2.16) 

as a function of N, for two values of g. The optimal values are seen to be iV* = 12 and N* = 5, 
in agreement with the estimate (2.10). 

Perhaps the most important problem in asymptotic analysis is how to go beyond optimal 
truncation, incorporating in a systematic way the small exponential effects (2.4). In this first 
section we will address this problem in the context of asymptotic series appearing in ODEs. 



2.2 Formal power series and trans-series in ODEs 

Very often, the solution to a physical or mathematical problem is given by an asymptotic series, 
and one has a systematic procedure to calculate this series at any desired order (like for example 
in perturbation theory in QM). The first question we want to ask is the following: can we 
calculate, at lest formally, non-perturbative effects of exponential type (i.e. like the one in (2.4)) 
which must be added to an asymptotic series? Since these terms are typically not taken into 
account in classical asymptotics, in order to include them we have to consider generalizations of 
asymptotic series. The resulting objects are called trans-series and they were first considered in 
a systematic way by Jean Ecalle in his work on "resurgent analysis" [53] . 

An important class of asymptotic series and trans-series are the formal solutions to ODEs 
with irregular singular points. The simplest example is probably Euler's equation 

= ± (2.17) 
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which has an irregular singular point at z = oo. There is a formal power-series solution to this 
equation of the form 

oo 

= E^T' a n = A- n nl (2.18) 

n=0 

This solution is an asymptotic series, with zero radius of convergence, since the coefficients grow 
factorially with n. It is easy to see that one can construct a family of formal solutions to this 
ODE based on y>o(z), 

<p{z) = <p (z) + Ce~ Az (2.19) 

where C is an arbitrary constant parametrizing the family of solutions. This is our first example 
of a trans-series. This example has three important properties: 

1. The term added in (2.19) is non-analytic at z = oo, and it goes beyond the standard 
solution in the form of an asymptotic series given by (2.18). 

2. The resulting formal expression has two small parameters, z and e~ A ^ z . 

3. There is a relation, already pointed out above in the context of optimal truncation, be- 
tween the "strenght" of the non-perturbative effect, given by A, and the divergence of the 
asymptotic series. Namely, A encodes as well the next-to-leading behavior of a n at large n. 

These properties are typical of formal trans-series and will reappear in many examples. 

A more complicated example of ODE with trans-series solution is the well-known Airy equa- 
tion, 

if" = xcp. (2.20) 

The solutions to this equation, called Airy functions, are ubiquitous in physics. We are now 
interested in formal power series solutions to this equation around x = oo. It is easy to see that 
one such a solution is given by 



oo 



^W = ^e-^/3^ v -3n/ 2 , (2 . 21) 

where 

Here, the coefficients grow as 

a n ~A- n n\, A = ~\- ( 2 - 23 ) 

What is the trans-series solution to this equation? As it is well-known, there is another, inde- 
pendent formal power series solution to the Airy equation, given by 

1 oo 

Z Bl (x) = -r^e 2 * 3/2 /3 Y^{-lTa n x~^\ (2.24) 

The general "trans-series solution" to the Airy differential equation is just a linear combination 
of these two formal asymptotic series, 

C 1 Z Ai (x) + C 2 Z Bi (x). (2.25) 
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Notice that Z\i(x) and Z-Qi(x) have different leading exponential behavior. Moreover, the relation 
between these behaviors is in accord with the third property noted above in the context of 
Euler's equation: the growth (2.23) of the coefficients suggests that the trans-series added to the 
asymptotic solution (2.21) should have a relative exponential weight of 

e 4 * 3/2 / 3 (2.26) 

as compared to (2.21), which is indeed the case. 

In the case of non-linear ODEs the structure of trans-series solutions is much richer: linear 
ODEs have trans-series with a finite number of terms, while in nonlinear ODEs they have an 
infinite number of terms. A class of important examples which are relevant in many physical 
applications are the Painleve transcendants. We will focus on the cases of Painleve I (PI) and 
Painleve II (PII). 

Example 2.1. Painleve I. The PI equation is 

u{k) 2 - -u"(k) = k. (2.27) 
6 

This equation appears in many contexts. In particular, it gives the all-genus solution to two- 
dimensional quantum gravity (see [49] for a review). There is a formal solution to this equation 
which goes like u(k) ~ yfn as k — > oo: 

»(»)(«) = ^t^n _ ^ _ _ ^ _ i|| K -¥ + ...). (2 . 28) 

n=0 v 7 

The trans-series solution to Painleve I is a one-parameter family of solutions to (2.27) which 
includes exponentially suppressed terms as n — > oo: 

oo oo 

„(«) = £cVO(k) = ^J2c e K~f e" MK5/4 e W(K), (2.29) 
where C is a parameter, the constant A has the value 

A = (2.30) 



and 

oo 

e W( K ) = J>, nK " 5 "/4 (2.31) 

?1=0 

are asymptotic series. Since we have introduced an arbitrary constant C in (2.29), we can 
normalize the solution such that = 1. We will refer to the above series u^(k) with £ > 1 
as the the ^-th instanton solution of PI, while the solution u^^k) will be referred to as the 
"perturbative" solution. 

It is easy to see that the -£-th instanton solutions in the trans-series satisfy linear ODEs. 
For example, one has for £ = 1, 

(k), (2.32) 

and from this one finds, 



5 s 75 b 341329 



15 



e ( 1 )( K ) = l + -^ z -'i — + . . . _ ( 2 _ 33 n 

64^/3 8192 23592960^/3 

Recursion relationships for the coefficients n^ >ri in (2.31) can be found in [62]. 
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Example 2.2. Painleve II. The Painleve II equation is 

u"{k) - 2u 3 (k) + 2ku(k) = (2.34) 

This equation is of fundamental importance in, for example, random matrix theory and non- 
critical string theory. It appears in the celebrated Tracy- Widom law governing the statistics of 
the largest eigenvalue in a Gaussian ensemble of random matrices [116], in the double-scaling 
limit of unitary matrix models [107] and of two-dimensional Yang-Mills theory [70], and it also 
governs the all-genus free energy of two-dimensional supergravity [85] . As in the case of PI, there 
is a formal solution to PII which goes like u(k) ~ y/n as k — > oo: 

m , , r- 1 73 10657 13912277 

« u « =vt w ir 1? > k-^oo. (2.35) 

16 K2 512 /si 8192 542888 kT 

One can consider as well exponentially suppressed corrections to this "perturbative" behavior 
and construct a formal trans-series solution with the structure, 

OO OO 

(«) = £cV*)(«) = ^^C £ K -¥e- MK3/2 e W(K), k — y co, (2.36) 



where 



and 



A = \ (2.37) 



e^( K ) = Y,u e ,nK- 3n/2 . (2.38) 

n=0 

As before, we normalize the solution with u±fi = 1- The perturbative part u^°\k) is given by 
(2.35). The instanton expansions can be easily found by plugging the trans-series ansatz in the 
Painleve II equation. One finds, for example, for the one-instanton solution, 

e (i)(«) = i _ l!«-3/2 + i^_ K -3 _ . . . . (2 . 39) 
v ' 96 18432 v 7 

After these examples, we can now give some more formal definitions (see [39]). Let 

<p' = f(x,<p), (2.40) 

be a rank n system of non-linear differential equations. We assume that f (x, <p) is analytic at 
(oo, 0). Let Aj, i = 1, ■ ■ ■ , n, be the eigenvalues of the linearization 

A = -(!^(oo,0)).. . (2.41) 

By using various changes of variables, one can always bring the system to the so-called normal 
or prepared form 

ip' = -A<p-^B<p + g(z,<p), (2.42) 

where 

A = diag(Ai,--- ,A n ), S = diag(^i,--- ,fi n ), (2.43) 
and by construction g(z, ip) = 0(|y?| 2 , z~ 2 ip). We also choose variables in such a way that Ai > 0. 
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If we write 



Example 2.3. Airy equation in prepared form. Define 

z = x 3/2 , (2.44) 

so that the Airy equation reads 

V"(z) = ^(z)-±<p>(z). (2.45) 

Let us consider the matrix 

*w - (r 3/ l 'A) ■ <™» 

V 1 Az 1 ~T Az/ 

u(z) = S^M*), = (^ } ) , (2.47) 

we find 

u'(z) = -Au-^Bu + g(z,u), (2.48) 

with 

The formal trans-series solution to (2.42) is of the form 

<p = ^ + ^ C k e- k - Az z- k ^ k , (2-50) 

keN™\{0} 

where 

C = (d,--- ,C n ) (2.51) 

are free parameters and 

C k = Cf 1 •••C*\ (2.52) 
The functions <£> an d V'k are formal power series, of the form 



Vk = Z>ki»*" n (2-53) 



n>0 

We will also denote 



|k| = J^jfei. (2.54) 

i 

Remark 2.4. For linear systems, like the Airy equation, all the trans-series <^ k vanish when 
|k| > 2, so the general trans-series solution is of the form 

n 

^ = (^ + ^Qe- A ^-ft^- (2-55) 
i=i 
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Figure 2: Three paths which lead to solutions of the Airy equation. 



2.3 Classical asymptotics and the Stokes phenomenon 

Stokes, by mathematical supersubtlety, transformed Airy's integrals... 

Lord Kelvin, "The scientific work of Georges Stokes." 

So far we have worked at a formal level and we have obtained formal solutions, in terms of 
asymptotic series, to ordinary differential equations. A natural questions is: what is the meaning 
of these formal power series? In order to answer this question, it is useful to consider in detail 
the case of the Airy equation and the Airy function, and to use an integral representation (we 
follow here the discussion in chapter 4 of [102]). Let us consider the integral 



(2.56) 



where 7 is a path which makes the integral convergent. Three such paths are shown in Fig. 2, 
but not all of them are independent, since 



7i + 72 + 73 = 0. 



(2.57) 



It is easy to see that the above integral (2.56) gives a solution to the Airy differential equation. 
We will now focus on the function defined by the path 71: 



Ai(x) 



2vri 



dze x 



(2.58) 



71 



which defines the Airy function Ai(cc). This function is analytic in the complex plane. We set 

x = re iK (2.59) 

and rescale the integrand 



ur 



1/2 



(2.60) 
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We find in this way 



1/2 r 

Ai(x) = — / e r 3/ V K «-« 3 /3) dn _ (2.6i) 

27T1 J yi 



We now want to study the behavior of this function for \x\ 3> 1, by using the saddle-point 
method. We then focus on the integral 

e AS K («) du ^2.62) 



71 



where 

There are two saddle points 
with "actions" 



u 3 



S K (u)=e 1K u--. (2.63) 
U R = e iK/2 , u L = -e iK/2 (2.64) 

S K (u RL ) = ±^exp(|^. (2.65) 

It is easy to see that the formal asymptotic power series around the saddle point u L is precisely 
(2.21), while the one around ti R is (2.24). We introduce the notation 

i? RL = Re (S K (u R > L )) , / RL = Im (S K (u R > L )) . (2.66) 

The paths of steepest descent (ascent) passing through these points are those where the function 
Re(S K (u)) decreases (respectively, increases) most rapidly, as we move away from the critical 
point. We will also denote by Jk' L the steepest descent paths through u RL , respectively. We 
show some paths of steepest descent and ascent in Fig. 3. 

Let us now study what happens as we change the angle k. For 

N < y (2-67) 

the path 71 can be deformed into a path of steepest descent through the saddle point at it L (for 
k = 0, 7r/3, the steepest descent paths are shown in Fig. 3.) We therefore have 

Ai(x) = / 7£ , \k\ <—, (2.68) 

which leads to the asymptotics, 

2vr 

Ai(x) - Z Ai (x), \k\ < —. (2.69) 

When 

\k\ = 2vr/3 (2.70) 

the steepest descent path coming from the saddle at u L runs right into the other saddle point. 
At this angle we have 

lm(S K (u L )) = lm(S K (u R )). (2.71) 
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_ 27T 57T 

t\ — — ur — 

3 6 

Figure 3: Saddle-point analysis of the integral (2.62) for different values of K. The red dot on the left 
is the critical point u L , while the black point on the right is the critical point m r . The continuous lines 
represent paths of steepest descent, while the dashed lines are paths of steepest ascent. 

Values of k for which this happens are called Stokes lines. This is the place where the second 
saddle might start contributing to the integral. In fact, for 

2vr 

— < \k\ < vr (2.72) 

O 

the contour 71 gets deformed into a steepest descent path passing through u together with a 
steepest descent path passing through u R , and 

Ai(x) = / 7£ +/ 7 «. (2.73) 

However, in this range of k the saddle u R gives an exponentially suppressed contribution to the 
asymptotics. In classical asymptotic analysis, subleading exponentials are not taken into account, 
and the asymptotics is still given by the contribution from u L : 

M(x) ~ Z Ai (x), \k\ < vr. (2.74) 

However, when x < 0, both saddles have the same real part 

Re(S K (u L )) = Re(S K (u R )). (2.75) 

A line where this occurs is called an anti-Stokes line. Therefore, both saddles contribute to the 
asymptotics, which is then given by a linear combination of the two trans-series and Z-q\ (the 
precise combination can be obtained by a more detailed analysis, see [102]). One finally obtains 
an oscillatory asymptotics: 

Ai(x) ~ 1 ' cos -|x| 3/2 - - , x<0, |rr|^oo. (2.76) 



vr V3 1 1 4 
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The fact that different asymptotic formulae hold on different directions for the same analytic 
function is called the Stokes phenomenon. From the point of view of saddle-point analysis, what 
is happening is that the saddle point which appeared on the Stokes line, at k = 2-7r/3, is no longer 
subdominant at k = tt, and it has to be included in the asymptotics. 




Figure 4: Saddle-point analysis of the Airy function Ai(x). Full lines (in red) are Stokes lines, while 
dashed lines (in blue) are anti-Stokes lines. On the Stokes lines k = ±27r/3, a second saddle appears 
in the integration contour. This saddle is subdominant when 2n/3 < \n\ < tt and does not contribute 
to classical asymptotics. However, at n = 7r, the saddle is not subdominant anymore and leads to an 
oscillatory asymptotics. 

The saddle-point analysis of the Airy integral is summarized in Fig. 4. There are Stokes 
lines at 

2-7T 

Stokes : arg(x) = 0, arg(:c) = ± — , (2.77) 

O 

and anti-Stokes lines at 

anti-Stokes : arg(x) = ir, =t — . (2.78) 

3 

In light of the example of the Airy function, we can now understand the meaning of the 
formal trans-series solutions to ODEs. ODEs have "true" solutions which are, generically, mero- 
morphic functions on the complex plane. The classical asymptotics of these solutions is given by 
particular trans-series solutions, i.e. by linear combinations of formal solutions to the ODE. Due 
to the Stokes phenomenon, this combination changes as we change the angular sector where we 
study the asymptotics. Therefore, formal trans-series solutions are the "building blocks" for the 
asymptotics of actual solutions. 

In the context of systems of ODEs, Stokes and anti-Stokes lines are defined as follows. Let 
us consider our system in prepared form (2.42). The directions where an exponential is purely 
oscillatory, i.e. 

Re(A i z) = (2.79) 

are called anti-Stokes lines. Along these directions, terms which used to be exponentially sup- 
pressed become of the same order than the leading term. This leads to an oscillatory asymptotic 
behavior. The directions where 

Im(A^) = 0, (2.80) 
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are called Stokes lines. These are the directions where subleading exponentials start contributing 
to the asymptotics. If we consider the Airy equation in prepared form, the eigenvalues are ±2/3. 
In terms of the variable z defined in (2.44), the Stokes lines are at arg(z) = 0,7r, while the anti- 
Stokes lines occur at arg(z) = ±vr/2, which, when translated to the variable x, give the structure 
found above. 



Figure 5: The sectors of the K-complex plane of opening 8ir/5 where the tritronquee solutions w_ 2 , 

uo, Mi, U2 are represented by the formal series u^(n). The dots in the remaining sector of opening 2n/5 

represent the infinite number of poles that the tritronquee solutions have there. 

Example 2.5. Painleve I and its tritronquee solutions. The meaning of the formal power series 
solution to PI (2.28) can be clarified by looking at actual solutions of the Painleve I ODE. It 
can be shown (see for example [60] and [81]) that there exist five different genuine meromorphic 
solutions of (2.27) with the asymptotic power expansion (2.28) in one of the sectors of the z- 
complex plane of opening 87r/5, see Fig. 5: 



Along the remaining sector of opening 2n/5, the asymptotics involves elliptic functions and 
the solutions have an infinite number of poles. 

2.4 Beyond classical asymptotics: Borel resummation 

So far we have discussed formal power series and trans-series. These formal power series give 
asymptotic approximations of well-defined functions which in the case of ODEs are their "true" 
solutions. Our next question is: to which extent can we recover the original, "non-perturbative" 
solution, from its asymptotic representation? Since asymptotic series are divergent, the answer 
to this question is not obvious. 

In fact, various answers have been proposed to this question. The more traditional answer 
is to use the optimal truncation procedure discussed in section 2.1. This gives a reasonable 
approximation to the original function in some regions of the complex plane, but it typically 
becomes a bad one in other regions. A nice illustration is provided again by the Airy function. 
Let us run the following numerical experiment, proposed by Berry in [21]. In optimal truncation, 
we approximate 



-14ji/5 



-631/5 




-2jt/5 



-23t/5 





(2.82) 
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where 



N* 



4/3|aj 



and [] denotes the integer part. Let us now plot the parametrized curve 

(Re(Ai(|x|e iK )),Im(Ai(|a;|e iK ))), < k < vr, 



(2.83) 



(2.84) 



in the complex plane, for fixed |x| = 1.7171, and let us compare it with the corresponding curve 
computed by using the r.h.s. of (2.82). The result is shown in Fig. 6. The approximation is quite 




Figure 6: Parametrized plot of the real and imaginary parts of the exact Airy function (red) compared 
to the optimal truncation of the asymptotic approximation (blue). The thin green line signals the angle 
K = 2tt/3. 



good in the region 



< arg(x) < 



2tt 



(2.85) 



and it becomes worst and worst as we approach arg(x) = n. The reason is simple: we are missing 
an exponentially small correction! This correction is born on the Stokes line arg(x) = ^ and 
becomes more and more important as we approach the anti-Stokes line, where it is of the same 
order than the term we are keeping. It is clear that, in order to reproduce the full function, 
we must find a way to incorporate the exponentially small corrections. Notice as well that, in 
optimal truncation, only a finite number of terms in the asymptotic expansion are actually used, 
and the remaining terms cannot be used to improve the estimate. 

The most powerful way to go beyond optimal truncation and solve the above problems is 
probably the technique of Borel resummation. Let 



ip{z) 



\ - Or, 



n>0 



be a factorially divergent series with a r , 



n\. Its Borel transfom is defined by 



£(C) = a n 



n>l 



_c_ 

(n-1)!' 



(2.86) 



(2.87) 
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This series defines typically a function which is analytic in a neighbourhood of the origin. If the 
resulting function can be analytically continued to a neigbourhood of the positive real axis, in 
such a way that the Laplace transform 

/•oo 

/ e-*0(C)d< ( 2 - 88 ) 
J o 

converges in some region of the z-plane, then the series ip(z) is said to be Borel summable in that 
region. In that case, 

poo 

s(ip)(z) = a + / e-*£(C)dC. (2.89) 
J o 

defines a function whose asymptotics coincides with the original, divergent series (f(z), and 
s(<p)(z) is called the Borel sum of f(z). 

Remark 2.6. Sometimes we want to perform the Borel resummation along an arbitrary direction 
in the complex plane, specified by an angle 6. It is then useful to introduce the generalized Borel 
resummation 

pe lS oo 

s e (<p)(z) = a + e-* c £(C)dC. (2.90) 

J o 

A crucial issue in the analysis of Borel resummation is the location of the singularities of 
y?(C)- It is easy to see that, if 

a n ~ A~ n n\, (2.91) 

then (p(C) is analytic in an open neighborhood of radius A around £ = 0. There is a singularity 
at z = A, which can be a pole or a branch point. If the singularity is not on the positive real 
axis, the integral (2.88) defining the Borel resummation is typically well defined and reconstructs 
the original function, see Fig. 7. 




Figure 7 : The Borel transform defines an analytic function in a neighbourhood of f = 0, of radius p = A. 
There is a singularity on the circle |£| = A, shown here as a red point on the negative real axis. If we can 
analytically continue this function to a neighbourhood of the positive real axis, and its Laplace transform 
exists, we say that the series is Borel summable. 

Example 2.7. Borel resummation and Euler's equation. Let us consider Euler's equation (2.17) 
with A = — 1. Then, the formal solution is the asymptotic series 

/ \ (-l) n n! 

= £ Sir-- ( 2 - 92 ) 

n>0 
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Its Borel transform is 

oo 1 

m = E(- 1 ) n " 1 c n_1 = r— r ( 2 - 93 ) 

n=l ^ 

Since ip(Q has no singularities on the positive real axis, we can define the Borel resummation 

8(<p){z) = P" 2 ^, (2.94) 

which defines an analytic function in the region 

Rez>0 (2.95) 
and reconstructs a true solution to the original differential equation. 



Example 2.8. Borel resummation and the c = 1 string. Let us consider the following asymptotic 
series, 

n>l v ' 



where B n are the Bernoulli numbers. Since l?2fc = for k > 1, this series has only even powers of 
z. This series appears often in string theory. It gives the genus expansion of the c = 1 string at 
the self-dual radius (see for example [84]), and it also appears in the asymptotic 1/JV expansion 
of the partition function of the Gaussian matrix model (see for example [92]). To see that the 
series is asymptotic, notice that 

B n = -2(-l) n / 2 (2vr)- n n!C(n), n > 2. (2.97) 

The zeta function behaves at infinity as, 

C(n) w 1, n -> oo, (2.98) 

up to exponentially small corrections in n, so the series is alternating and factorially divergent. 
Its Borel transform can be computed explicitly, 



(2.99) 



It has singularities along the imaginary axis, at the points C = 27rmi, m = ±1, ±2, • • • . It has no 
singularities along the positive real axis, and the Borel transform 



(2.100) 



gives the Borel resummation of the original series for z > (see [107] for more details on this 
and related examples). 

Example 2.9. Borel resummation and the Airy function. In the case of the Airy function we 
can proceed as follows. Let us define 



vai(*) = E5. (2.101) 



z 

n=0 
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where the coefficients a n are given in (2.22). Its Borel transform can be explicitly computed as 
a hypergeometric function 

au(0 = ~ (2 ' 102) 

and it has a branch point singularity at £ = —4/3. The Borel resummation, 

POO 

s (ip M ) (2) = a + / (?Ai(C)e~ 2C dC (2.103) 

JO 

is well-defined if 

Re(z) > (2.104) 

and it reconstructs the full Airy function in the region (2.104) (see for example [45]). We then 
have, 

Ai(x) = —L^e-^/'s (ip Ai ) (z), x = z 2 '\ (2.105) 

2x i / 4 V7T 

In terms of arg(x), this representation is valid as long as 

|argx| < J. (2.106) 



Figure 8: The paths C± avoiding the singularities of the Borel transform from above (respectively, below). 

Very often one encounters asymptotic series whose Borel transform has singularities on the 
positive real axis. In this case one needs some prescription in the integral (2.88) to avoid the 
singularities. A standard procedure to do this is to consider lateral Borel resummations. Let Cg± 
be a path going from to +e ie oo and avoiding the singularities of (p(Q on the direction with 
angle 9 from above (resp. below). For 9 = 0, we will simply write the paths as C±, and in this 
case they have the form shown in Fig. 8. The lateral Borel resummations are then defined as, 

se± fa) (z) = a + ( e -<£(C)dC, (2.107) 

JC e ± 

provided the integral is convergent. Notice that, even if the original series has real coefficients 
and we choose a direction along the real axis, since the lateral Borel resummations are computed 
by integrals along paths in the complex plane, they lead in general to complex-valued functions. 
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Example 2.10. Lateral Borel resummations and Euler's equation. Let us consider the Borel 
transform (2.93) along the negative real axis. Since there is a singularity at £ = —1, we are 
forced to perform lateral Borel resummations: 

s„±(<p)(z)= [ e-<-^-. (2.108) 

These integrals give two different solutions to the original differential equation for z < 0. Their 
difference can be computed as a residue, 

S7r+ {if) (z) - *„_ (if) (z) = f f^dC = 27riRes f= _ 1 f^ : = 2vrie 2 , (2.109) 

and it is exponentially small along the negative real axis. In fact, it is the exponentially small 
term appearing in the trans-series solution (2.19). See [110] for a detailed discussion of this 
example. 

Example 2.11. Lateral Borel resummations for the Airy function. Let us consider the Borel 
transform (2.102) along the negative real axis in the z variable. Notice that, in terms of the x 
variable, this corresponds to arg(x) = 27r/3, i.e. to a Stokes line. Notice that, along this direction, 
the coefficients of the series defining the Airy function (2.21) are no longer alternating, and this 
is the standard indication that the series is not Borel summable along such a direction. Indeed, 
since there is a singularity at ( = — 4/3, we are forced to perform lateral Borel resummations: 

s„± (VAi) (z) = a + ( £ Ai (CK 2C dC. (2.110) 

JCn± 

It can be shown that (see, for example, [45]) 

Sn+ (¥>Ai) (z) - (<p M ) (z) = -ie 42/3 s^ (0 Bi ) (z) (2.111) 

where 

oo 

^BiW = E(- 1 ) n 5 = ^i(-^). (2.H2) 

n=0 

Notice that, since (fBi(z) is analytic on the negative real axis, the above Borel resummation is 
well-defined. We find that, as in the previous example, the difference of lateral resummations is 
given by a trans-series solution. As we will see in a moment, this is the way in which the Stokes 
phenomenon manifests itself in the context of Borel resummations. 

Let us now see in general how to apply Borel resummation to the study of ODEs. If yo(z) is 
a formal solution to an ODE, its Borel resummations (or lateral Borel resummations) will provide 
functions which solve the ODE and have the asymptotic behavior given by yo(z). However, we 
can add exponentially small corrections to this solution without changing the asymptotics. In 
general, the multi-parameter family 

<P±(z; C±) = a± (cp ) (z) + £ C^- k ^e- A - k2 S± (cp k ) (z), (2.113) 

k 

which is obtained by doing lateral Borel resummations on the formal trans-series solution (2.50), 
is a good solution for sufficiently large \z\ which asymptotes to yo(z), provided the non- vanishing 
terms in the trans-series are such that 

Re(A-kz)>0. (2.114) 
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The reciprocal is true: any solution to the ODE can be represented by such a Borel-resummed 
trans-series for an appropriate choice of C's. This is one of the main consequences of Ecalle's 
theory of resurgence, see for example [110, 32, 39] for detailed statements and proofs. We see that 
the main advantage of the Borel resummed version of asymptotic analysis is that we can make 
sense of small exponentials, i.e. we can incorporate the information encoded in the trans-series 
in a systematic way. This is not the case in classical asymptotics. 

What is the interpretation in the context of Borel resummation of the Stokes phenomenon? 
In classical asymptotics, Stokes lines indicate the appearance of small exponentials, as we have 
seen in the analysis of the Airy function: we pass from (2.68) to (2.73). However, this jump 
is only noticed in the classical theory when we reach the anti-Stokes line. Once we use Borel 
resummation, we can give a "post-classical" version of the Stokes phenomenon. Let us consider 
a Stokes direction, which we take for simplicity to be arg(z) = 0, corresponding to the eigenvalue 
A = Ai > 0. Along this direction, there are two families of solutions y±(z;C±), obtained by 
lateral summations from below and from above. By uniqueness we should expect these two 
solutions to be related. Indeed, we have the relation 

¥>+(*; C) = ¥>_(*; C + S), (2.115) 

where 

S = (5i,0,--- ,0) (2.116) 

is called the Stokes parameter associated to the Stokes line arg(z) = 0, and it is an imaginary 
number when the coefficients of the trans-series are real. At leading order in the exponentially 
small parameter exp(— Az) we find 

¥>(>;+(*) " ¥>(>;-(*) « Si^^V(ia.,o),(4 (2-117) 

The relation (2.115) is the Borel-resummed version of the Stokes phenomenon. It says that the 
coefficients of the trans-series solutions have discontinuous jumps along the Stokes direction. The 
results (2.109) and (2.111) are two particular examples of this general result, where the Stokes 
line is the negative real axis for the z variable, and the Stokes parameters in these examples are 
S = 2iri for the Euler equation, and S = — i for the Airy function. For a pedagogical explanation 
of (2.115) in the case of ODEs, see [110, 32]. The generalization to systems of ODEs is presented 
in [39]. 

The relationship (2.117) has a very nice interpretation in terms of Borel transforms, which 
we will make explicit for simplicity in the case of a first order ODE with /3 = 0. Let us denote 
the perturbative solution by (fo(z), with the form (2.86), and the first trans-series by 



n>0 



The l.h.s. of (2.117) can be written as 

'dCe-^(C), (2.119) 



where 7 is homotopic to the contour C+ — C_, and encircles the singularities of the Borel transform. 
Then, (2.117) tells us that the structure of <p(C) around the singularity at ( = A is of the form 

^{A + = Si (g| + + holomorphic) . (2.120) 
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where ^i(C) is the Borel transform of (2.118). This is easy to check: the integral (2.119) can be 
evaluated by using (2.120). The first term in (2.120) gives the residue at the pole ( = A, namely 

,Sie~ A Vi,o, (2.121) 

while the second term in (2.120) gives the integral of the discontinuity of the log, namely 

POD _ 

J ° n>l Z 

We then reconstruct the l.h.s. of (2.117). 

The relationship (2.117) is then telling us that the singular behavior of the Borel transform 
of the perturbative series is related to the first instanton trans-series. This shows that "per- 
turbative" and "non-perturbative" phenomena are intimately related, at least in this example. 
Later on we will see that this relationship has a powerful corollary, namely it gives an asymp- 
totic formula for the large order behavior of the coefficients of the perturbative series. Equation 
(2.120) is an example of the resurgence relations discovered by Jean Ecalle, and it forms the basis 
of the so-called "alien calculus" of his theory (see [31] for an introduction). In fact, (2.120) is 
just the tip of the iceberg, since relations of this type connect all the formal power series in the 
trans-series, and not only the perturbative series and the first instanton. 

Example 2.12. In the case of the Airy function, one can use the relation (2.111) to derive the 
following formula for the Airy function along the negative real axis [45] 

Ai(x) = ^={e- 2 ^ /2 / 3 ^ /2 fcAi) (z) +ie 2 * 3/2 / 3 ^ /2 ( m ) (*)}, M*) < 0, (2-123) 

where x = z 2 / 3 . This is the exact version of the oscillatory asymptotics given in (2.76) and (2.73). 

Example 2.13. Painleve II and the Hastings-McLeod solution. As an example of how to re- 
construct a "true" function from the asymptotics in the case of non-linear ODEs, including 
exponentially small corrections, we consider the example of Painleve II, whose formal structure 
was discussed in Example 2.2. This equation has a Stokes line at arg(fc) = 0. The PII equation 
has a solution, called the Hastings-McLeod solution, «hm(k) which is uniquely characterized by 
the following properties: 

1. It is real for real k. 

2. As k —7- co it asymptotes 

«hm(«) ~ k 1/2 - (2.124) 

3. As k — >■ -co it asymptotes 

u(k) ~ e - 2 ^(- K ) 3/2 /3. (2.125) 

This solution of PII plays a crucial role in the Tracy- Widom law [116] and in the double- 
scaling limit of unitary matrix models [42]. The Hastings-McLeod solution is shown in Fig. 9. 
Since it is a solution to Painleve II, at sufficiently large k one must be able to express it as the Borel 
resummation of the trans-series (2. 36). The Borel transform has a singularity at ( = A = 4/3, 
therefore we have to use lateral resummations along the positive real axis. The relation (2.115) 
reads in this case, 

u + {k-C) = U-(k;C + S) (2.126) 
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-4-2 2 4 

Figure 9: The Hastings-McLeod solution to Painleve II. As k — > oo it asymptotes s/k, while as n — > — oo 
it decays exponentially as indicated in (2.125). 

where the subscripts ± refer to the two lateral resummations of the full trans-series solution 
(2.36), and 




(2.127) 



is the Stokes parameter. By using results on the non-linear Stokes phenomenon for the Painleve 
II equation [60] one can show that [94] 

uhm(k) = u+(k;-5/2). (2.128) 

Notice that this solution is real. Complex conjugation changes the sign of 5 and also exchanges 
the integrals along the contours C + and C_ , so we have 

u + (k; -5/2)* = u_(k; 5/2) = u+(k; -5/2) (2.129) 

where we used (2.126). The connection to Borel resummed formal solutions gives the correct 
"semiclassical" content of the Hastings-McLeod solution, and one easily shows that 

u + (k; -5/2) = hv® + uL 0) ) - ^5 2 (nf + u^) + ■■■ (2.130) 
2 o 

The equation (2.115) is very important conceptually, when interpreted from a physical point 
of view. As we will show in detail in these lectures, the formal power series yo has often the 
interpretation of a "perturbative" series, while the trans-series yk have the interpretation of 
non- perturb at ive effects. The coefficients C give the strength of these effects. But one clear 
implication of (2.115) is that this strenght is not well-defined unless we give a prescription to 
perform the Borel resummation of the power series appearing in the formal solution. Indeed, 
different prescriptions lead to different coefficients, and in particular the upper and lower lateral 
resummations differ by a shift involving the Stokes parameters. Notice that there is a "compen- 
sation" effect, in the sense that we can change simultaneously the resummation prescription and 
the strenght of the non-perturbative effects so that the solution is unchanged. This is indeed 
the content of (2.115). This phenomenon was noticed in the literature on renormalons in gauge 
theories (see for example [73] for a clear presentation), as well as in the study of instantons in 
QM [125], and is called in those contexts the cancellation of non-perturbative ambiguities. 
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2.5 Non-perturbative effects and large order behavior 

An important consequence of the Borel resummation technique is a relationship between the 
asymptotic behavior of the coefficients of a perturbative series, and the first instanton or trans- 
series solution. This type of relationships were anticipated in [50] and in the work on the large 
order behavior of quantum perturbation theory [90], and it will be important in the following 
lectures. We have already seen in the examples of the Euler equation and the Airy function that 
the "action" appearing in the trans-series controls the large order behavior of the perturbative 
coefficients. This is a general phenomenon. We will now give a heuristic derivation of an asymp- 
totic formula for the coefficients, in the case of a first order ODE with (3\ = (see [41, 65] for 
details and rigorous proofs). We will also write down various generalizations of the result. 
We will denote the "perturbative" series and its Borel transform by 

Mz) = T,%&> ft(0 = E^C". (2.131) 

Therefore 



Z — TV 

n>0 n>0 



^±1 = J_ / dC Mg. (2.132) 

where Co is a contour around the origin. We know that tfo{() has a singularity at £ = A, and 
we can deform the contour Co so as to enclose this singularity (there are singularities at other 
points with |£| > A, but they give exponentially small corrections as compared to what we are 
computing). We can then write 

^±1 ~ ± [ dC?#. (2.133) 

We know the singularity structure of f>o(C) near £ = A thanks to the result (2.120): there is a 
pole with residue —Si <pifl/2m, and a logarithmic discontinuity. Changing variables £ = A + £, 
we obtain 

a n +i Sm,o / 1 dC Si ^ <p ltk 1'°° C^ 1 i a r91 o,v 



n 



2vri J ((A + C) n+1 2vri 2vri ^ (fc - 1)! 7 (^ + C)™ +1 

t -,-,rwr(„ t i- t ) 



Since 

/ (A + C)"+ lUS "" r(n + l) 

we have the following result for the all-order asymptotics of the coefficients: 

■■-^^ n-'f. (2 - 136) 

This is a beautiful result. It implies that the leading asymptotics of the original ("perturba- 
tive") series is encoded in the first trans-series (or "one-instanton" ) solution (there are further, 
exponentially suppressed contributions associated to the higher singularities). Conversely, all 
the information about the formal one-instanton series is encoded in the asymptotics of the per- 
turbative series. Notice that the leading and next-to-leading order of the asymptotics is given 
by 

a n ~ n\A- fl (2.137) 
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as in the examples discussed above. However, (2.136) contains much more information. In 
particular, the Stokes parameter Si plays a crucial role in the asymptotics, and in fact (2.136) 
provides a method to determine numerically this parameter in cases in which it is not known 
analytically. 

Before considering generalizations of the above equation, let us work out in some detail an 
interesting example which will illustrate most of the considerations of this first lecture. 

Example 2.14. A Riccati equation. Following [110, 26], let us consider the ODE 

^ = ¥ ,-i(6- + 6V)- (2-138) 

dz z 

Here, fr^ are real constans, and we assume that 

(3 2 = -b-b + (2.139) 

is positive. The ODE (2.138) generalizes the Euler equation (2.17), which is obtained (for A = 
— 1) when b + = 0. Equation (2.138) is a particular case of the so-called Riccati equation, which 
is characterized by being quadratic in the unknown function. It is easy to see that there is a 
formal solution of (2.138) around z = oo which is given by 

¥*>(*) = 6 "E!£' ( 2 - 14 °) 

n>l 

where the coefficients a k are obtained from the non-linear recursion 

fe-i 

a k+ i = -ka k - P 2 ^2 a e a k ^ e , a = 1. (2.141) 

l=\ 

Explicitly, we find 

The Riccati ODE (2.138) has a full trans-series solution of the form 

v(z,C) = Y,C n e nz cp n (z), (2-143) 

n>0 

therefore there is a series of "multi-instantons" with action nA, and A = —1. The Stokes 
parameter for this ODE has been computed exactly in [26], and it is given by 

Si = -2vri6_a(/3), a(p) = ^M. ( 2 .144) 

Tip 

Using (2.136), we obtain the following asymptotics for the coefficients a n : 

a n ~ a((3)(-l) n -\n - 1)! (l + O (±\ V (2.145) 

This is easy to test numerically. To do that, one simply studies the asymptotic behaviour of the 
sequence 
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which should converge towards the constant value cr(/3). However, the resulting convergence is 
quite slow. One can accelerate it by using Richardson transformations. Let us assume that a 
sequence s n has the asymptotics 



oo 



(2-147) 



n 

k=0 



for n large. Its iV-th Richardson transform can be defined recursively by 



,(o) 



(TV) _ (AT-1) n ( (N-l) _ (N-l), N > 1 

The effect of this transformation is to remove subleading tails in (2.147), and 



(2.148) 



1 

n 1 ' 



^-co + ^r^r). (2.149) 



The values s^ give numerical approximations to cq, and these approximations become better 
as N, n increase. Once a numerical approximation to Co has been obtained, the value of c\ can 
be estimated by considering the sequence n(s n — Co), and so on. In Fig. 10 we plot the original 
sequence s n and its first and second Richardson transforms, for (3 = 1/3. The convergence 
towards 

a « 0.826993... (2.150) 

(2) 

is quite fast, and s^q gives an approximate value for this constant which agrees with the right 
value up to the seventh decimal digit. 




The asymptotic result (2.136) has many generalizations. For example, when f3 ^ 0, and the 
trans-series have the structure 

oo 

W (z) = z-? k Y,^nZ- n , (2.151) 

n=0 
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a simple generalization of the above argument shows that 

an „ |lr(n - P)A~ n+ P Y . ^ kAk . (2.152) 

One can also generalize the argument to higher order ODEs. In that case, as seen in (2.113), 
there are various possible instanton actions for the trans-series. The asymptotics is governed by 
the trans-series which correspond to the smallest actions, in absolute value. If there are instanton 
actions which have the same smallest absolute value but different phases, the asymptotics is 
obtained by adding their different contributions, see [41] for a precise mathematical statement. 
For example, for the coefficients u n $ of the perturbative solution (2.28) of PI, one has the 
following asymptotics [81, 79] 

Mo „~ A-^lT(2n - \) Ml + V -= ^ 1, (2.153) 

which comes from two instanton actions ±A. In this expression, u± i are the coefficients of the 
1-instanton series €^(k) appearing in (2.31), A is the instanton action (2.30), and 

S, = -i A (2.154) 

is a Stokes parameter. 

Since the trans-series solutions are also formal, asymptotic series, one can ask what is the 
asymptotic behavior of their coefficients. In the same way that the asymptotics of the pertur- 
bative coefficients is encoded in the first trans-series, it turns out that the asymptotics of the 
coefficients of a given trans-series is encoded in higher trans-series solutions. The study of this 
question requires the full machinery of resurgent analysis, see [65, 62, 10] for various results along 
this direction. 

2.6 Lessons 

We can now summarize some of the results that can be learned from the study of asymptotic 
series appearing in ODEs. All of these results will have a counterpart when we look at instanton 
corrections in quantum theories and string theories, and therefore they constitute a sort of "role 
model" for their study. 

1. The standard perturbative contribution (yo(x), in the context of ODEs) is a factorially 
divergent, asymptotic series. At the formal level, one can also obtain trans-series solutions. 
These will correspond to perturbation theory around instanton solutions, i.e. to non- 
perturbative effects. 

2. The weights of the instanton solutions are a priori undetermined. Therefore, the general 
trans-series solution gives a multi-parameter family of formal solutions. This is the non- 
perturbative ambiguity. 

3. By Borel resummation, the family of formal solutions becomes a family of "true" solutions. 
Therefore, we obtain a family of non-perturbative completions. If we have a non-perturbative 
definition of the theory we can fix the non-perturbative ambiguity by choosing the values 
of the parameters that reproduce the non-perturbative definition. Along directions where 
the series is Borel summable, the original solution is typically reconstructed by Borel re- 
summation of the perturbative series solely. 
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4. Along Stokes lines there are different prescriptions for resummation, due to singularities 
in the Borel transform. Their difference is purely non-perturbative and defines the Stokes 
parameter. This is the "resurgent" version of the Stokes phenomenon. The reconstruc- 
tion of the non-perturbative solution involves in a crucial way the Borel-resummed non- 
perturbative effects, as we showed in Example 2.13 for the Hastings-McLeod solution of 
PII . 

5. The large order behavior of the perturbative expansion along a Stokes line encodes the 
action of the instanton, the Stokes parameter, and the coefficients of the first instanton 
correction. This relation is extremely powerful, since it says that the large order behavior 
of perturbation theory knows about non-perturbative corrections. In cases where there is 
no clear technique (or even framework!) to address the computation of non-perturbative 
effects, the large order behavior of the perturbative series gives an important hint about 
their structure. 



3. Non-perturbative effects in Quantum Mechanics and Quantum Field Theory 
3.1 Trans-series in Quantum Mechanics 

Before discussing non-perturbative effects in QFT and matrix models, it is instructive to first 
consider simple quantum-mechanical examples, where the analysis of non-perturbative effects 
can be made in detail. We will focus on the ground state energy of one-dimensional particles 
with Hamiltonian 

H = £+V(q). (3.1) 
We will assume that the potential V(q) is of the form 

V{q)= l -q 2 + gV int {q), (3.2) 

where gV\ n t(q) is the interaction term. A typical example is the quantum anharmonic oscillator, 
where 

V{q) = ^ + \x\ (3.3) 

The ground state energy of a quantum mechanical system in a potential V(q) can be com- 
puted in a variety of ways. The most elementary method is of course Rayleigh-Schrodinger 
perturbation theory, and the resulting series has the form 



E = J2»ng\ (3.4) 



n=0 

where (in units where h = 1) 

a = - (3.5) 

is the ground state energy of the harmonic oscillator, and we have an infinite series of corrections 
due to the interaction term in (3.2). In order to make contact with QFT it is instructive to 
calculate the series (3.4) in terms of diagrams. To do this, we first notice that the ground state 
energy can be extracted from the small temperature behavior of the thermal partition function, 

Z(0) =tre-^ (/3) , (3.6) 
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S(q) = / dt 

J -aii 



\{m? + V{q{t)) 



(3.9) 



as 

E = - lim \ (3.7) 
In the path integral formulation one has 

Z(/3) = j V[q(t)]e- S ^, (3.8) 

where S(q) is the action of the Euclidean theory, 

'■a/2 
-a/2 

and the path integral is over periodic trajectories 

q(-(3/2) = gr(jS/2). (3.10) 

The path integral defining Z can be computed in standard Feynman perturbation theory by 
expanding in Vi Q t(q). We will actually work in the limit in which {3 — > oo, since in this limit 
many features are simpler, like for example the form of the propagator. In this limit, the free 
energy will be given by (3 times a /3-independent constant, as follows from (3.7). In order to 
extract the ground state energy we have to take into account the following facts: 

1. Since we have to consider F(/3) = log Z(j3), only connected bubble diagrams contribute. 

2. The standard Feynman rules in position space will lead to n integrations, where n is the 
number of vertices in the diagram. One of these integrations just gives as an overall factor 
the "volume" of spacetime, i.e. the factor (3 that we just mentioned. Therefore, in order to 
extract E(g) we can just perform n — 1 integrations over R. 



t' 




9 
4 



Figure 11: Feynman rules for the quantum mechanical quartic oscillator 
For (3 — > oo the propagator of this one-dimensional field theory is simply 

dp e [ P T e"l r l 



(3.11) 
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For a theory with a quartic interaction (i.e. the anharmonic quartic oscillator) 

Mnt(g) = \q A (3.12) 

the Feynman rules are illustrated in Fig. 11. One can use these rules to compute the perturbation 
series of the ground energy of the quartic oscillator (see Appendix B of [17] for some additional 
details). We have, schematically, 

a n = ^2 (connected vacuum bubbles) . (3.13) 

For example, the diagrams contributing up to order g 3 are shown in Fig. 12, and after performing 
the integrals over the propagators one finds, 



1 3m_2im J , 333 

2 4V4/ 8 V4 



which agrees with the result of standard Rayleigh-Schrodinger perturbation theory. 

A basic property of the above perturbative series is that the coefficients a n grow factorially 
when n is large. Moreover, this behavior is due to the factorial growth in the number of diagrams 
(the Feynman integrals over products of propagators only grow exponentially). To see this, 
remember that a n can be computed as a sum over connected quartic graphs. The total number 
of connected graphs with n quartic vertices is given by 

i((x 4 )«)W, (3.15) 

where 

f°° dxe~ x / 2 x 4n 

{{X T) = £ dse-*V2 (3 - 16) 

is the Gaussian average. By Wick's theorem, the average counts all possible pairings among n 
four-vertices. The superscript (c) means that we take the connected part of the average. Since 



we find 

-<(* 4 r> - (4n " 1)!! - (4re)! (3 is) 

n\ i[X)) - n\ "4"n!(2n)!- [6A *> 

As n — > oo this behaves like 

4 2n n!, (3.19) 

i.e. there is a factorial growth in the number of disconnected diagrams. One could think that 
there might be a substantial reduction in this number when we consider connected diagrams, 
but a careful analysis [16] shows that this is not the case: at large n, the quotient of the number 
of connected and disconnected diagrams differs from 1 only in 0(1/ n) corrections. We conclude 
that there are ~ n! diagrams that contribute to a n . The resulting factorial behavior of the 
perturbative series of the quartic oscillator can be verified by a detailed consideration of Feynman 
diagrams [19] (see [15] for a review of these early developments). Therefore, we conclude that 
the perturbative series for the ground state energy is a formal, divergent power series. This series 
gives at best an asymptotic expansion of the true non-perturbative ground-state energy, defined 
in terms of the exact spectrum of the Schrodinger operator. 
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Figure 12: Feynman diagrams contributing to the ground state energy of the quartic oscillator up to 
order g 3 . 

We can now ask what is the analogue of the trans-series for this type of problems. As it is well- 
known (see for example the discussion in the textbook [124]), there are instanton contributions 
to the thermal partition function. For concreteness, let us consider again the quartic oscillator 
and let us suppose that the coupling constant is negative, i.e. g = —A, with A > 0, so that we 
have a potential of the form shown in the left hand side of Fig. 13. In this case, the Euclidean 
action has non-trivial saddle-points. The EOM reads 

-q(t) + q(t)-Xq 3 {t) =0. (3.20) 

In the limit j3 — > oo one finds the following trajectory with E = 0, 

^ = <^ h ^hor (3 - 21) 

where, for /3 finite, 

-P/2 <t < /3/2 (3.22) 

is a free parameter. When (3 — > oo, such a trajectory starts at the origin in the infinite past, 
reaches the zero of the potential V(q) at t = to, and returns to the origin in the infinite future, 
i.e. 

q c — > 0, t — > ±oo, 



t = to, 




(3.23) 



As is well-known, the Euclidean action can be regarded as an action in Lagrangian mechanics 
with an "inverted" potential —V{q), and the non-trivial saddle-point described above is simply 
a trajectory of zero energy in this inverted potential. 
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Figure 13: The inverted potential relevant for instanton calculus in the quartic case. The instanton or 
bounce configuration q c (t) leaves the origin at t = — oo, reaches the zero (2/A)s at t — t , and comes back 
to the origin at t = oo. 

The partition function around this non-trivial saddle can be computed, at one loop and 
for large (3, by using standard techniques, which we will not review in here (a very complete 
and updated discussion can be found in chapter 39 of [124]). It can be seen that this saddle- 
point is unstable: it has one, and exactly one, negative mode. This means that the instanton 
contribution is imaginary. A detailed analysis, which can be found in for example [38, 124], 
shows that this one-instanton calculation determines the discontinuity of the partition function 
for negative values of the coupling: 



discZ(-A) = Z(-\ + ie)- Z(-X-ie) = 2iImZ(-A) 



(3.24) 



This leads to a discontinuity in the ground-state energy, as a function of the coupling. In the 
case of the quartic oscillator one finds, at one loop, 



discE(-X) = 2iIm#(-A) 



8i 



^/2ttA 



-A/A 



where 



.4 



4 
3 



(3.25) 



(3.26) 



is the action of the saddle (3.21) for A = 1. This imaginary correction to the energy has a clear 
physical interpretation: since for negative coupling the potential is unstable, a particle in its 
ground state will eventually tunnel. The width of the ground state energy 



r = 2|ImJ5|, 



(3.27) 



is inversely proportional to the life-time of the ground state. 

The above calculation is just the one-loop approximation to the one-instanton sector. But if 
we consider multi-instanton expansions at all loops, we expect to find for the ground state energy 
a trans-series structure of the form 



where E^ is the asymptotic, perturbative series (3.4), and 



oc e 



-lA/g 



(3.28) 



(3.29) 
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are the ^-instanton corrections, themselves asymptotic expansions in the coupling constant z. In 
the case of the quartic oscillator they have the structure 



= //V M /*X>nA (3-30) 

n=0 

which is identical to (2.151) (here the expansion is around z = 0, while in (2.151) we do the 
expansion around z = oo). In other cases, like the double- well potential analyzed in [123, 124, 125] 
they are more complicated and include terms of the form z n \ogz. 

The structure of the trans-series in QM suggests that the instanton action determines the 
positions of the singularities of the Borel transform, and that the large order behavior of the 
coefficients in the perturbative series (3.4) is controlled by the first instanton contribution. These 
expectations are indeed true, and one can show that much of the structure appearing in ODEs can 
be extended to the analysis of quantum-mechanical potentials in one dimension. In particular, 
the "resurgent" structure of the formal trans-series calculating the energies of bound states in 
QM has been established, following the work of Voros [117], in [46]. 

In the case of the quartic oscillator, the Borel transform of the formal power series (3.4) 
has a singularity at £ = —4/3. It is therefore Borel-summable along the positive real axis (i.e. 
for g > 0), and its Borel resummation is indeed the exact ground-state energy, as defined by 
the Schrodinger operator. This was originally proved in [69], and a proof using the theory of 
resurgence can be found in [46]. To understand the large order behavior of the series (3.4) we 
have to take into account the presence of the instanton at negative g = —A < 0. In this case, 
one has to consider lateral resummations of E(g) along the negative real axis. The discontinuity 
disc.E(— A) gives then the difference between lateral Borel resummations, and the result (3.25) 
can be interpreted as the analogue of (2.117) in the theory of ODEs: the asymptotic expansion 
of this difference is given (at leading order) by the first instanton correction to the energy, which 
has the general structure 

s+(E)(-z) - 8-{E)(-z) « SiE^(z), (3.31) 
where Si is a Stokes parameter. At one- loop we have the asymptotic result, 

s+(E)(-z) - s-(E)(-z) « — =e~^. (3.32) 



Notice, in particular, that the coefficient of the one-loop calculation of the instanton partition 
function gives the Stokes parameter of the problem. 

As in the theory of ODEs, one can use this result to derive the large order behavior of the 
coefficients a n in (3.4). Writing 



SiE^ (-z) = iz^- A ' z c n z n , (3.33) 

n=0 

we find the asymptotic growth 

27T { 1=1 U m =i(n ~ P ~ m) J 
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This formula is exactly like the one in (2.152), with the only difference of the extra factor (— 
which is due to the fact that we do the perturbative expansion in the variable g = —z. Plugging 
in the concrete values of the quartic oscillator for the different quantities, i.e. 

P = ~\, c ° = 4 \/!' ^ = 4 /3- (3.35) 
we find for the large-order behavior 

y/6 /3 



The story of this result is a fascinating chapter of modern mathematical physics (see [114]). 
The behavior of the a n at large n was first obtained by Bender and Wu in [17] by studying 
numerically the sequence of the first seventy-five coefficients. They were even able to guess, from 
these numerical experiments, the exact form of the prefactor in (3.36). In a subsequent (and 
classic) paper [18], they showed that the result (3.36) could be derived analytically by looking at 
the one-instanton sector. 

In the case of the quartic oscillator with positive coupling g, the origin leads to a stable 
quantum-mechanical ground state. The instanton solution only appears when one inverts the 
sign of the coupling, and as a result the perturbative series for positive g is alternating and Borel 
summable. In other cases, like for example the cubic oscillator 

V{x) = l -x 2 - gx 3 , (3.37) 

the origin is always unstable quantum-mechanically for any real value of the coupling constant 
g. An elementary calculation shows that the action of the instanton mediating the tunneling is 



^ A=* 
g^ 15 



2 , A=-. (3.38) 



The large order behavior of the coefficients in the formal power series for the ground state energy 

oo 

E(g) = Y,a n g 2n (3.39) 

n=0 

can be computed by using a small modification of (3.34) -essentially, one has to consider the 
discontinuity along the positive real axis of g 2 and then there is no sign alternating factor. 
Therefore, 

a n ~ A~ n V (n + 1/2) (3.40) 

where A is given in (3.38), see for example [5] for a derivation of this result. In this case the 
series is not Borel summable, reflecting the instability of the perturbative ground state. 

In general, in one-dimensional quantum-mechanical problems, we will have complex instanton 
solutions with complex actions. They lead to perturbative series which are Borel summable and 
have an oscillatory character. As we mentioned in the case of ODEs, the large order behavior is 
controlled by the instantons with the smallest action in absolute value, and the phase of the action 
determines the oscillation period of the series. Let us analyze in some detail a very instructive 
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example, following [29] (a useful discussion can be also found in chapter 42 of [124]). Let us 
consider a particle situated at the origin of the potential 



9 



V{gx), V(x) = l -x 2 - jx 3 + ^x 4 . 



The ground state energy has the expansion 

E {g) = ^2 an9 



2n 



n>0 



In this example there are two different situations (see Fig. 14): 

1. For I7I > 1, the origin is not an absolute minimum, which is in fact at 



XQ 



3j+ J-8 + 9Y 



2. For I7I < 1, the origin is the absolute minimum. 



(3.41) 



(3.42) 



(3.43) 
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Figure 14: On the left: the potential V(x) in (3.41) for 7 = 0.95, where the origin is the absolute 
minimum. On the right: the potential (3.41) for 7 = 1.01; the origin is now unstable to quantum 
tunneling. 

In the first case |7| > 1, the vacuum located at the origin is quantum-mechanically unstable, 
and there is a real instanton given by a trajectory from x = to the turning point 



x + = 7 - 

The action of this instanton can be written as 
A r*+ 



Sr. 



r 



A = 2 I dx(2V(x)) 
lo 



r 



--+7 2 -- 7 (7 2 -l)log^. 
3 ' 2 ,w ; 7-1 



The one-loop prefactor for this instanton, appearing in (3.33), is given by 



7T 



(3.44) 
(3.45) 

(3.46) 



The behavior when I7I < 1 is obtained by analytic continuation of this instanton configura- 
tion, which is now complex. In fact, there are two complex conjugate instantons described by a 
particle which goes from x = to 

x = j±iy/l- 7 2 . (3.47) 
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We have then to add the contributions of both instantons. Since cq becomes imaginary when 
|7| < 1, adding the complex-conjugate contributions of the two instantons gives 

a k ~ T(k + \/2)lmA' k - 1 / 2 . (3.48) 

More generally, if we have a quantum-mechanical problem involving a complex instanton and its 
complex conjugate, and 

A=\A\e~ ie \ Co = | Co |e i0c , (3.49) 
the large order behavior, obtained by adding the contribution of the two instantons, is oscillatory 

a k ~ r(fc - (3)\A\~ k+ P cos ((A; - (3)6 A + C ) . (3.50) 

As in the case of ODEs analyzed in section 2, when a perturbative series is Borel summable 
(like in the case of the quartic oscillator with g > or in the potentials with complex instantons) , 
the Borel resummation of the perturbative series reconstructs the non-perturbative answer. There 
are two types of situations where there is no Borel summability: the first case corresponds to 
perturbative series around unstable minima, like the quartic oscillator with g < or the cubic 
oscillator. A different situation occurs in the case of the double-well potential. In that case, 
there is a stable ground state but the perturbative series is not Borel summable, and one has 
to consider lateral Borel resummations. The ground state energy can be reconstructed from the 
Borel-resummed perturbative series and the Borel-resummed instanton or trans-series solutions, 
in a way which is similar to the analysis of the Hastings-McLeod solution of Painleve II in 
Example 2.13, see [123, 46, 125] for more details on this quantum- mechanical problem. 

3.2 Non-perturbative effects in Chern Simons theory 

The extension of the above results to realistic QFTs is plagued with serious difficulties. Probably, 
the most important one is the fact that in renormalizable QFTs there are other sources of 
factorial divergence in perturbative series, namely renormalons (see [20]). Renormalons are 
particular types of diagrams which diverge factorially due to the integration over momenta in 
the Feynman integral. Due to the existence of renormalons, the analysis of the large order 
behavior of perturbation theory inspired by QM does not extend straightforwardly to standard 
QFTs. There are however QFTs where renormalon effects are absent, like Chern-Simons (CS) 
theory and many supersymmetric QFTs, and we will focus here on this "toy" QFTs, and more 
particularly on CS theory, where many different aspects of non-perturbative effects are relatively 
well understood. 

CS theory is a QFT defined by the action 

S = —?- I TrfAAdA+^-AAAAA). (3.51) 

Here, A is a G-connection on the three-manifold M, where G is a gauge group. We will mostly 
consider G = U(N), and in this case our conventions are such that A is a Hermitian N x N 
matrix-valued one-form. Gauge invariance of the action requires [47] 

keZ. (3.52) 

The 3d QFT defined by this action is a remarkable one: it is exactly solvable, yet highly nontrivial, 
and provides a QFT interpretation of quantum invariants of knots and three-manifolds [121]. 
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The partition function of the theory on M is defined by the path integral 

Z{M) = j VAe iS ^ (3.53) 

and in principle it can be computed by standard perturbative techniques (see [96] for a review). 
The saddle-points of the CS action are just flat connections 

F(A) = 0. (3.54) 

These are in one-to-one correspondence with embeddings 

m(M)^G, (3.55) 

where G is the gauge group. In principle, the path integral (3.53) has various contributions 
coming from expansions around the different saddle-points. The perturbative sector is defined 
by expanding around the trivial flat connection 

A = 0, (3.56) 

while instanton sectors are associated to non-trivial flat connections. Formal expansions around 
these instanton sectors define the analogue of trans-series for this QFT. 

Example 3.1. Lens spaces. The lens space L(p, 1) has fundamental group ni(L(p, 1)) = Z p . 
The set of U(N) flat connections is given by homomorphisms 

Z p -> U(N), (3.57) 

modulo gauge transformations. These are in turn given by splittings of U(N) into p factors 

U(N) -> £/(iVi) x U(N 2 ) x • • • x U(N P ), (3.58) 

corresponding to the homomorphism 



£ -> diag | • • • , f-y- ,e~\ 



(3.59) 



where 

{ = exp . (3.60) 

Therefore, instanton sectors are in one-to-one correspondence with partitions of iV into p nonzero 
integers. The CS action evaluated at the flat connection labelled by {Nj}j=\ t ... jP is given by 

A= 7 f±( j -lfN j . (3.61) 

We can now ask what is the nature of the perturbative series appearing in CS theory. It 
turns out that, generically, perturbation theory around the trivial connection A = is factorially 
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Figure 15: AS and IHX relations. 

divergent. The reason for this is the same as in QM, namely, the factorial growth in the number 
of diagrams. Let us see this in some detail. We will denote by 



F(M,g,g s 



(3.62) 



the contribution of the trivial connection to the free energy log Z of CS theory. Here, g is the 
Lie algebra associated to G, and 

2vri 



9s 



k 



(3.63) 



Using standard perturbative techniques, it is easy to see that the free energy can be written as 
a formal power series of the form 



F(M,g, 5s ) = ^ £ c r (M)W s (T)g^. 



(3.64) 



n=1 re^£ c) (0) 



Let us spell out in detail the ingredients in this formula. We first construct the space of Feynman 
diagrams, ^( c )(0). This is the space of connected, trivalent diagrams with no external legs (i.e. 
connected vacuum bubbles) modulo the so-called IHX and AS relations, shown in Fig. 15. It is 
graded by the degree of the diagram n, which is half the number of vertices (and also equals the 
number of loops minus one) : 

^ C )(0) = 0- =1 4 C) W- (3-65) 

A basic fact is that, for each n, this space has finite dimension. The very first dimensions are 
listed in Table 1. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


d{n) 


1 


1 


1 


2 


2 


3 


4 


5 


6 


8 



Table 1: Dimensions d(n) of An \$) up to n = 10. 
An explicit choice of basis up to n = 5 is shown below: 
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n = 2: 




(3.66) 

The second ingredient is the weight system. This is an instruction to produce a number for 
each diagram, given the data of a Lie algebra with structure constants and a Killing form, 

[T a ,T b ] = f abc . (3.67) 

To each trivalent vertex we associate the structure constant f a b c as shown in Fig. 16. In QFT 
we call this "computing the group factor of the diagram IY' The final ingredient is cp(M). It is 
simply given by the Feynman integral associated to the graph T. It is possible to show that each 
cr(M) is a topological invariant of M. For example, 

cg(M) = Casson invariant of M. (3.68) 




Figure 16: Weight system. 

Based on our experience with the quantum anharmonic oscillator, we should ask how many 
diagrams we have at each loop order n + 1. It has been shown by Garoufalidis and Le in [63] 
that 

dim (^ c) (0)) ~ n!, n > 1. (3.69) 

Therefore, the series (3.64) will be factorially divergent. Interestingly, there is no other source of 
factorial divergences. These divergences could come from the weight factors, or from the Feynman 
integrals. However, it is easy to see that the weight factors can only grow exponentially. It is 
also shown in [63] that the Feynman integrals grow with the degree as 

|cr(M)|~C^, (3.70) 
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where Cm is a constant that depends on the three-manifold under consideration. In QFT terms 
this means that Feynman integrals grow at most exponentially, i.e. that there are no renormalons 
(since in a renormalon diagram with n loops, the Feynman integral diverges itself factorially, as 
nl). 

Example 3.2. Chern-Simons theory on Seifert spheres. The general prediction of factorial 
divergence can be verified in detail for CS theory on Seifert homology spheres. A Seifert homology 
sphere is specified by r pairs of coprime integers (pj, qj), j = 1, • • • , r, and is denoted by 

M = x(^,--- ,^). (3.71) 



One also defines 



r 



p = Up» H = p Ef- ( 3 - 72 ) 

1 1 Ps 

s=l s=l 

H is the order of the first homology group Hi(M, Z). When r = 1, Seifert spaces are just lens 
spaces, 

L(p,q)=X(q/p). (3.73) 

The partition function of CS theory on a Seifert space can be written as a matrix integral. This 
was shown for G = SU{2) by Lawrence and Rozansky [88], and it was extended to arbitrary 
gauge groups in [91]. There are two types of non-trivial flat connections: the reducible ones, 
and the irreducible ones. In a Seifert space and when G is a simply-laced group, reducible flat 
connections are labelled by elements in 

t e Ar/HA r , (3.74) 

where A r is the root lattice. The contribution of such a connection to the partition function is 
given by (up to an overall normalization, see [91] for the details) 

Z(M) oc / dAe -AV2a.- fc t-A ni=in tt >o2sinhfe > (3 75) 

J rL>o(2smh^) 

where 

P 

9s = jj9s, (3.76) 

A belongs to the weight lattice A^, a > are the positive roots, and the products are computed 
with the standard Cartan-Killing form. In the case of U(N) we have rather 

t G Z N /HZ N , (3.77) 

and we write 

N N 
i=l i=l 

where is the orthonormal basis of the weight lattice, and < U < H — 1. We then find, 

Z(M) oc / dA, e~W s ££i S£i ^ IE=i n,<, 2 sinh ^ ^ 
' *=i II, . ; (2sinh A 2 A ' 
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The integration contour in (3.75), (3.79) is chosen in such a way that the Gaussian integral 
converges. 

The case of SU(2) is particularly simple. Up to an overall constant, the contribution of the 
trivial connection to the partition function is just an integral, 



Z<x Jd\e- x2 /^f(X;p s ), f(X; Pj ) = ^sinh fj ^2 sinh . (3.80) 
This can be expanded in power series in g s , 



n=0 

and one finds [88] 



f( 2n )(0) 

Z^y^a n g:, a n = J - ^, (3.81) 



oo 



p \ n 

<*, (3-82) 
i.e. we have a factorial divergence, as expected. The growth is controlled by 

A = -^f. (3.83) 

We would expect this quantity to be the action of a non-trivial saddle-point of the theory, and 
indeed this is the action of an irreducible flat connection on the Seifert manifold. 

The results of [63] about the structure of the CS perturbative series rely on a mathematical 
construction for this series called the LMO invariant [89], which has been much studied during 
the last years. The structure of the general trans-series, which correspond to the perturbative 
expansion around a non-trivial instanton solution, is less understood, and no mathematical con- 
struction has been proposed so far. Questions on classical asymptotics and Borel summability 
in CS theory have started to be addressed only recently, see [61, 122] for some results and/or 
conjectures. 

3.3 The l/N expansion 

The problem of non-perturbative effects and asymptotics becomes much more interesting when 
we look at gauge theories in the l/N expansion [115]. In this expansion, the free energy and 
correlation functions of the gauge theory are expanded in powers of N or of the coupling constant 
g Sl but keeping the 't Hooft parameter 

t = g s N (3.84) 

fixed. For example, the expansion of the free energy around the trivial connection (i.e. what we 
have called the perturbative series) is re-organized as 

oo 

F = Y J F g {t)9l 9 -\ (3-85) 

9=0 

where F g (t) is a sum over double-line graphs or fatgraphs of genus g. In this reorganization of 
the theory, the dominant contribution comes from the genus zero or planar diagrams. 

In the case of CS theory, the structure of the l/N expansion can be made very explicit, as 
follows. Consider a graph T in .4,^(0), and apply the thickening rules depicted in Fig. 17 and 
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Figure 17: Thickening an edge. 




Figure 18: Thickening a marked vertex. 

Fig. 18. The thickening rules can be regarded as a map that associates to each diagram T a 
formal linear combination of fatgraphs T g ^, which are Riemann surfaces with boundaries and are 
classified topologically by their genus g and number of boundaries h: 

r->£^( r ) r M- (3-86) 

9,h 

It is easy to see that the weight system of U (N) can be written very explicitly in terms of 
fatgraphs [43, 13], 

W <N) (T) = ^p g<h (T)N h . (3-87) 

9,h 

An example is shown in Fig. 19. One then finds the following expression for the free energy 
around the trivial connection: 

oo 

F(M,n(N),g s ) =J2J2 c ^ M )P9A r ) Nh 9f (r) - V{r \ (3-88) 

9=0 T gth 

where E(T), V(T) are the number of edges and vertices in T (these topological data do not depend 
on the fattening of the graph). If we now use Euler's relation, 

E(T) -V(T) = 2g -2 + h (3.89) 

we see that F(M, u(N),g s ) is given by the formal series (3.85), where F g (t) is defined as a formal 
infinite sum over all fatgraphs Tg^ with fixed g 

F g(t) = J2 a 9M\ a 9 ,h = 5>r(M)p ff , ft (r). (3.90) 

As is well-known (see for example the classic review in [37]) the 1/N expansion described 
above can be implemented in any U(N) gauge theory where the fields transform in the adjoint 
representation of U (N) and it can be applied to any gauge-invariant observable of the theory. 
The structure of the free energy as a double power series, see (3.85) and (3.90), which we have 



e- 

Figure 19: Fatgraphs obtained from the theta diagram. 

written above based in the analysis of CS theory, can be easily seen to hold in any theory with 
fields in the adjoint of U(N). The 1/N expansion is particularly clean in theories where the 
coupling constant g s does not run, i.e. in conformal field theories and topological field theories. 

The first question that we have to ask in the search for non-perturbative effects is: what 
is the nature of the formal power series appearing in the theory, like for example in the series 
defining the free energy of the theory expanded around the trivial connection? In the case of the 
1/N expansion, since there are two parameters, we have two different questions to ask: 

1. What is the nature of the formal power series in t appearing in (3.90), defining F g (t)l 

2. For a fixed t, what is the nature of the power series in g s appearing in (3.85)? 

The quick answer to these questions is the following. In theories with no renormalons, the 
functions F g (t) are analytic at the origin, i.e. the power series (3.90) have a finite radius of 
convergence which, moreover, is common to all g. However, for fixed t, the functions F g (t) grow 
like 

F g (t) ~ (2g)\ (3.91) 

and the series (3.85) diverges factorially. 

The analyticity of the 1/N expansion at fixed genus was first noticed in [87] and analyzed 
in some models by 't Hooft [115]. It can be proved in detail in simple U(N) gauge theories, such 
as matrix models (see [66] for a recent study) and CS theory [64]. The factorial growth of the 
1/N expansion was pointed out, in a slightly different context, in [112]. 





Figure 20: Counting flowers with n petals (figures courtesy of [48]). 

The basic reason for the analyticity of the F g (t) is that the number of fatgraphs of fixed 
genus grows only exponentially, and not factorially. A nice example of this contrasting behaviour 
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is the counting of flower-like graphs. The number of flowers made out of n "thin petals", like the 
one shown on the left in Fig. 20, is given by the number of possible contractions in a vertex with 
2k legs, and equals (2n — 1)!!, which grows factorially. In contrast, the number of "genus zero" 
flowers with n "thick petals," like the one shown on the right in Fig. 20, is given by the Catalan 
number (see for example [48]) 

C n = , ^ , ~4" (3.92) 
[n + iy.nl 

which only grows exponentially. 

We will now given an argument for the analyticity of the free energy in theories without 
renormalons, following [64]. Let us come back to the formal power series appearing in (3.90). It 
is easy to see that 

P g ,hF) < C 2 p 9 - 2+h . (3.93) 

For example, in a theory with a purely cubic interaction, like CS theory, each vertex gives two 
resolutions (see Fig. 18), and the maximum number of terms is 2 V . Since in a theory with cubic 
interactions we have 

3V = 2E (3.94) 

we deduce 

P g ,h( r ) < 2V = ^ E ~ V = ^ 9 ~ 2+h - (3.95) 

In QCD there are also quartic vertices, which from the point of view of this counting can be 
regarded as two cubic vertices joined by an edge, leading to a similar estimate. The next step is 
to analyze the Feynman integrals, cr- If the theory has renormalons, they can grow factorially 
with the number of vertices. But in a theory without renormalons they grow only exponentially 
in the number of vertices, and we can write 

|cr| ~ C 2 F 9 ~ 2+h . (3.96) 

This has been shown to be the case for a large class of diagrams in QM [15], and it has been 
proved to be the case in CS theory, by using the formulation in terms of the LMO invariant [63]. 
We then have, 

a g , h ~ (C p C F ) 29 - 2+h N g , h , (3.97) 

where N g ^ is the number of double-line diagrams with genus g and h holes, counted with the 
appropriate weight. For example, if we normalize all vertices of degree p with a factor l/p\, the 
weight of a diagram F g h is given by 

1 

|Aut(F M )|' (3 - 98) 

i.e. the inverse of the order of its automorphism group, see [48]. The counting of fatgraphs 
(weighted by their automorphism group, as above) has been developed very much both in com- 
binatorics and in mathematical physics. The main result we have in this respect is that 

Ngfi-CZC^pgy., (3.99) 

see for example [64]. We conclude that 

a g , h ~ (2g)\C{Cl (3.100) 

Therefore, for fixed genus, and assuming the growth (3.96), the functions F g (t) defined by the 
power series (3.90) are analytic at t = with a finite radius of convergence p common to all g. 
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This is the analiticity result we wanted to establish. Generically p < oo, and there is typically a 
singularity t c in the t-plane somewhere in the circle of radius p, 

\t c \=P- (3.101) 

This argument shows as well that for fixed t (and inside the common domain of convergence) the 
sequence F g (t) will diverge like (2g)\. This was our second claim above. 

The analyticity of the genus g free energies F g (t) suggests that t should be regarded naturally 
as an arbitrary, complex variable. In the same way, g s should be also thought of as a complex 
variable, although in the original gauge theory it can only take special values: in Yang-Mills 
theory g s is real and positive, while in CS theory it is of the form (3.63). The complexification 
of these variables is natural in the context of some large N dualities. For example, in dualities 
between large N CS theory and topological string theory [68], the 't Hooft parameter is identified 
with a complexified Kahler parameter of a Calabi-Yau manifold. We will then take the point of 
view that all gauge theory parameters belong to a complex moduli space. 

3.4 Large N instantons 

We have seen that the power series defining free energies F g (t) at fixed genus are analytic at t = 0. 
Therefore, there is in principle no indication of non-perturbative effects to be taken into account: 
if one can perform an analytic continuation of the function F g (t) to a region in the complex plane 
which includes all physical values of the 't Hooft coupling, then one has solved the problem. This 
procedure is sometimes called weak-strong coupling interpolation, since one starts with a function 
defined at weak 't Hooft coupling (i.e. in a neighbourhood of t = 0) and ends up with a function 
defined for some angular region in the complex plane where \t\ 3> 1. This analytic continuation 
to the physical region can be very often performed in interesting QFTs, and it is unique (we will 
see below an example of this, concerning the free energy of ABJM theory). 

However, as we have also seen, the price to pay for the analyticity of F g (t) for fixed g is that, 
at fixed t, the sequence of free energies grows doubly-factorially with the genus. More precisely, 
one has the growth 

F g {t) ~ (2g)\(A(t))- 2 ^ 5 »1, (3.102) 

where A(t) is a function of t. This behavior has been found in many actual, simple models of the 
1/N expansion, like matrix models (and their double-scaling limits) and in some supersymmetric 
gauge theories. 

As we have seen, in theories without renormalons, the standard factorial growth of the 
perturbative expansion is typically related to the existence of instantons in the theory. It is 
then natural to suspect that the growth (3.102) of the 1/N expansion should be also due to 
instanton-like objects which we will call large N instantons, and that A(t) is the action of such 
an object. 

What is a large N instanton? In general, large N instantons are built upon classical instan- 
tons. To see this, consider a gauge theory with coupling constant g s and 't Hooft parameter t, 
as well as an instanton solution whose action (including the coupling) is given by 

S c = —. (3.103) 

9s 

We will also assume that A is of order one at large N: 



A~C?(1). 



(3.104) 



For example, in Yang-Mills theory the usual instantons with low instanton number satisfy this 
property. This is due to the fact we can build an instanton by using just an SU(2) subgroup 
of U(N). This is also the case in matrix models, where, as we will see in detail, instantons are 
obtained by eigenvalue tunneling of one eigenvalue out of N. There are instanton configurations 
whose action is of order N ("giant instantons"), but we will not consider them here. 

We can now do perturbation theory around such an instanton configuration. The one-loop 
fluctuations give a term with the generic form (at large N) 



(3.105) 



where c\N is the number of zero modes, or collective coordinates of the instanton, at large N. 
This factor comes from the canonical normalization of the modes in the path integral, since we 
can always normalize the fields in such a way that the action has an overall power oil/g s . Putting 
both things together and expressing everything in terms of g s and the 't Hooft parameter, we 
find 



c\N 



-A/g s 



exp 



9s 



where 



A(t) = A-dtlog - )+0(t) 



(3.106) 



(3.107) 



is called the large N instanton action. It is given by a series in t which incorporates, on top 
of the classical action and the one-loop fluctuations which we have written down explicitly, all 
vacuum, connected planar diagrams (at all loops) in the background of the classical instanton 
action. To see how these appear, let us focus for simplicity on the interaction given by a cubic 
vertex, as in Fig. 18. Let us consider fluctuations around the instanton solution A 



A = A + A'. 



The action for the fluctuations will include a vertex of the form 



(3.108) 





Figure 21: The instanton vertex (3.109) (left) and a planar diagram contributing to the large N instanton 
action A(t) a term of order t 3 (right). 



EWj(4(V< ( 3 - 109 ) 

and involving the instanton background. We can represent this vertex in the double-line notation 
as in Fig. 21, where the red line ending on the blob corresponds to the instanton background. 
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It gives a factor of g s , but only the interior line gives a factor of N after tracing over. A simple 
example of a diagram contributing to the instanton action is the one depicted on the r.h.s. of 
Fig. 21. The inner closed lines gives a factor of iV 3 , and the diagram is proportional to 

Tr(A 3 )N 3 g 2 s = -t 3 Tr(A'), (3.110) 

9s 

since there are nine edges E = 9 and seven vertices V = 7, so the power of gf~ v is two. This 
diagram gives a correction of order i 3 to A(t). We see that the calculation of large N instanton 
actions in realistic theories is of course difficult, since we have to sum up an infinite number of 
planar diagrams (in the same way that calculating the genus zero free energy involves adding up 
an infinite number of diagrams at all loops). 

An alternative, more general way to think about large N instantons is in terms of large N 
effective actions. It is believed that gauge theories at large ./V can be reformulated in terms of 
a "large N effective action" with coupling constant (or h constant) equal to 1/N, and involving 
a field sometimes called "master field" [120]. In this effective theory, correlation functions at 
large ./V are obtained simply by solving the classical equations of motion of the effective action 
in the presence of sources. A large ./V instanton is simply an instanton solution of this large 
N effective theory, i.e. a saddle point with finite action. This is in general different from the 
usual instanton configurations, which are saddle point of the classical action. However, as we 
have seen, large N instantons can often be thought of as deformations of the classical instantons, 
where the deformation parameter is the 't Hooft parameter: as it is manifest in (3.107), when 
t — > we recover the action of the "classical" gauge theory instanton. Explicit examples of large 
N instantons were obtained in the CP^ model in [1, 103], as deformations of classical instantons. 
A particularly beautiful example is the large N instanton of two-dimensional Yang-Mills theory 
obtained in [70]. In the next section we will see an example of a large N instanton in a simple 
toy model, namely matrix quantum mechanics. 

Large N instantons play an important role in the structure of the 1/N expansion. First of 
all, in the calculation of physical quantities in a large N theory, we expect to have a trans-series 
structure with a perturbative sector around the trivial instanton sector, and then a series of 
sectors corresponding to 1/N expansions around large N instantons. The instanton sectors are 
weighted by 

e- A( * )/Ss . (3.111) 

As it will be obvious in examples, A(t), the large N instanton action, is in general a non-trivial 
function of the 't Hooft parameter t. As expected from the examples developed in these lectures, 
this is the same quantity appearing in (3.102) and controlling the subleading large g asymptotics. 

The second important property of large N instantons is their dynamical role in triggering 
phase transitions. If Re(A(t) / g s ) > 0, large N instantons are suppressed exponentially at large 
N (or small g s ). This might lead to think that "instantons are suppressed at large N" but as 
Neuberger pointed out in [104], this is not necessarily the case. It might happen for example that 
A(t) vanishes at a particular value of t, and in this case the contribution of instantons become as 
important as the perturbative contributions. The value of the 't Hooft parameter for which A(t) 
vanishes signals very often a large N phase transition, or a critical point, in the theory. This 
behavior is nothing but the large N version of the jump in the asymptotics occurring along an 
anti-Stokes line. It turns out that the critical value of the 't Hooft parameter is also, in many 
cases, the first singularity t c in the t-plane which we found in (3.101). One of the first examples 
of such a transition was found in [72, 118], and it was argued in [104] that this should be due to 
the vanishing of the instanton action, see [94] for a detailed verification. 
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Remark 3.3. Large N instantons as D-branes. Notice that the diagrams contributing to the 
large iV instanton action are similar to the diagrams involving an external boundary, coming from 
a Wilson loop for example. This means that, if the large N theory has a string dual, the large N 
instanton should be associated to a hole in the worldsheet -in other words, to a D-brane. The 
relation between large N instantons and D-branes in string theory goes back to the observation 
by Shenker [112] that the large order behavior (3.102) should be typical of genus expansions 
in string theory [112]. This observation was based on the duality between doubly-scaled large 
N matrix models (which also display this (2g)\ growth) and noncritical strings. According to 
Shenker, however, this growth should be a universal feature of string perturbation theory. The 
non-perturbative effects associated to this stringy behavior were later on identified with D-brane 
and membrane effects in [108, 14]. 

Example 3.4. A large N instanton in CS theory. Let us consider U{N) CS theory on the 
lens space L(2, 1), and the instanton (i.e. the non-trivial flat connection) characterized by the 
splitting 

(TV -1,1). (3.112) 

The path integral comes with a weight e~ A ^ s , where g s = 2g s (see (3.76)), and one finds from 
(3.61) that 

A=y- ( 3 - 113 ) 

Using the matrix model representation of the partition function given in (3.79), it is possible to 
calculate the large N instanton action built on this classical instanton [97]. The result is A(t)/g s , 
where 

A(t) = 2Li 2 (e" i / 2 ) - 2Li 2 (-e~*/ 2 ) = ^ - tlog (jj + 0(t 3 ), (3.114) 

which has indeed the expected structure (3.107). 

As we see, given the complexity of obtaining large N instanton actions in closed form, it is 
hard to make progress on the detailed understanding of non-perturbative effects in general large 
iV theories. In the next section we will make a detailed study on large N instantons in a simple 
setting, namely matrix models, but we will now present another instructive solvable example, 
namely large N instantons in Matrix Quantum Mechanics. 

3.5 Large N instantons in Matrix Quantum Mechanics 

In order to illustrate the general considerations on \/N expansions and large N instantons 
explained above, we will now discuss Matrix Quantum Mechanics (MQM), introduced and first 
studied in [28] . This is a quantum- mechanical model where the degrees of freedom are the entries 
of a Hermitian N x N matrix M and the Euclidean Lagrangian is given by 



L M = TV 



M 2 + V(M) , (3.115) 



.2 

where V(M) is a polynomial in M. Notice that this problem has a U(N) symmetry 

M UMU^ (3.116) 

where U is a constant unitary matrix, and it promotes the standard one-dimensional QM problem 
to a problem where fields are in the adjoint representation of a U(N) symmetry group and can 
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then be studied in the 1/N expansion. We will assume that the potential V(M) is of the form 



V(M) = iM 2 + Mnt(M) 



(3.117) 



where Vi n t(M) is the interaction term. The Feynman rules are the same as in the case of QM, 
with the only difference that we will now have "group factors" due to the fact that M is matrix 
valued. 
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Figure 22: Feynman rules for matrix quantum mechanics. 
The propagator of MQM is 



-M 



and for a theory with a quartic interaction 

Mnt(M) 



4iV 



:M 4 



(3.118) 



(3.119) 



the Feynman rules are illustrated in Fig. 22. The factor of N in (3.119) is introduced in order to 
have a standard large N limit, as we will see in more detail later. 

One can use these rules to compute the perturbation series of the ground state energy of 
MQM, which is obtained by considering connected bubble diagrams, as in conventional QM. Each 
Feynman diagram leads to a group factor which depends on N, i.e. each conventional Feynman 
diagram gives various fatgraphs that can be classified according to their topology. A fatgraph 
with V vertices and h boundaries will have a factor 

t v N h-v = t v N 2-2g^ (3 120 ) 

since the number of edges is twice the number of vertices, E = 2V and 

h + E-V = h-V. (3.121) 
Planar diagrams, as usual, are proportional to iV 2 , and the ground state energy has the structure 



E(t,N) = Y,N 2 - 2 °£ g (t). 



(3.122) 
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The first few terms in the expansion of £o(t) can be easily computed in perturbation theory, 

£ ( t ) = - + -t-—t 2 + —t 3 + --- (3.123) 
uw 2 8 256 1024 v J 

As first found in [28], the planar ground state energy in MQM can be obtained exactly by 
using a free fermion formulation. This exact result resums in closed form all the planar diagrams 
of MQM contributing to the ground state energy. This goes as follows. After quantization of the 
system we obtain a Hamiltonian operator 



H = Tr 

where 



i d 



2 



2dM 2 



+ V(M) , (3.124) 



ab 

In order to study the spectrum of this Hamiltonian, it is useful to change variables 

M = UMJ\ (3.126) 

where 

A = diag(Ai,A 2 ,--- ,Ajv) (3.127) 

is a diagonal matrix. It can be shown that, when acting on singlet states (i.e., states that are 
invariant under the full U(N) group), the differential operator (3.125) has the form, 



o=l a b^a 



— Tr 



2 9M 2 2 9A2 1 A 6 - A a dA a 



Due to the residual Weyl symmetry, a singlet state is represented by a symmetric function of the 
N eigenvalues, 

tt(Ai). (3.129) 
We now introduce a completely antisymmetric wavefunction 

$(A) = A(A)*(A), (3.130) 

where 

A(A) = H(X a - X b ) (3.131) 

a<b 

is the Vandermonde determinant. It is now easy to see that the original problem of calculating 
the energies for singlet states becomes the problem of calculating the energy of N non-interacting 
fermions (since the function (3.130) is completely antisymmetric) in an external potential V(X). 
We will now assume that this potential V(A) has good large scaling properties. More precisely, 
we will assume that the A-dependence of the potential V(A) is such that 

V(\) = Nv(-^=) , (3.132) 



where v(X) does not contain N. For example, the quartic potential considered above, 



- 50 - 



has the right scaling properties. This can be interpreted as saying that 

t = Ng s (3.134) 



is the 't Hooft parameter of the model, which is kept fixed as N — > oo. After rescaling A — > y/NX, 
we find that the one-body fermion problem reduces to 

tr <h +v(\)\<t> n {X) = e n <f> n (\) (3.135) 



2N 2 dX 2 



where 

e„ = —E n (3.136) 

and E n are the energy levels in the original one-body problem. The ground state energy is given 
by 

N n 

E(t, N) = Y,E n = NY,e n = N 2 £ (t) + ■■■ (3.137) 

n=l n=l 

To calculate So, we note that the quantum effects in (3.135) are controlled by h/N. Therefore, 
large N is equivalent to H small, and in the large ./V limit we can use the semiclassical or WKB 
approximation. In particular, we can use the Bohr-Sommerfeld formula to find the energy 
spectrum at leading order in h/N. We will write this semiclassical quantization condition as 

NJ(e n ) = n ~\, n>l, (3.138) 

where 

1 fMe) 



J(e) = — / dA v / 2(e- v(X)) (3.139) 

4(e) 

and Ai ) 2(e) are the turning points of the potential. If we denote 



?= n ^, (3-140) 



we see that (3.138) defines implicitly a function e(£) through 

J(e(0)=e (3.141) 

At large N, the spectrum becomes denser and denser, the discrete variable £ becomes a continuous 
one £ G [0, 1], and the sum in (3.137) becomes an integral through the rule 

N 



n=l J ° 



One then finds, 

So = I d£e(0- (3.143) 
J o 

We also define the Fermi energy of the system by the condition 

J(e F ) = 1. (3.144) 



- 51 - 



After some simple manipulations, one finds the expression 

»A 2 (e F ) 



So = e F 



— I 

3irh Jx 



dX 



\i(e F ) 



2(e F - V(X)) 



(3.145) 



In the case of the quartic potential, the ground state energy can be explicitly computed as a 
function of the 't Hooft parameter. We first obtain the Fermi energy, which is defined by (3.144). 
For the quartic potential, the integral J(e) can be easily computed in terms of elliptic functions. 
Let us denote, 

2 _V4^TT-i 



a 



yJAet + 1 + 1 



t ' t 
The turning points of the potential are ±a. We also introduce the elliptic modulus 

„2 



k 2 = 



a 2 + b 2 ' 



We then find 



1 



J(e) = -^(2t)5( a 2 + b 2 )2 b 2 K(k) + (a 2 - b 2 )E(k) 
dirn 



(3.146) 



(3.147) 



(3.148) 



where E(k), K(k) are complete elliptic integrals. The condition (3.144) defines the Fermi energy 
e F (t) as an implicit function of the 't Hooft parameter. The planar free energy is given by 



£ (t) = e F (t)-—l(t,e F (t)) 



(3.149) 



and involves the integral 



I(t,e) 



du 



( a 2_ n 2 )(6 2 + u 2 



2 

35 



(3.150) 



Va 2 + 6 2 {2(a 2 - 6 2 )(a 4 + 6a 2 b 2 + b 4 )E(k) + b 2 {2b 4 + 9a 2 6 2 - a 4 )K(/fc)}. 
The final result can be easily expanded in powers of t, and one finds 

So(t) 



17t 2 75t 3 3563t 4 

+ tt^tt - -^T^r + 0(t 5 ). 



1 t 

2 + 8 256 ' 1024 32678 



(3.151) 



The first few terms are in perfect agreement with the calculation in planar perturbation theory 
(3.123). 

One important remark on this result is that £o(t) is an analytic function of t at t = 0, in 
accordance with the general result explained above. This follows from the explicit expression 
for £q in terms of elliptic functions. The radius of convergence of the expansion (3.151) can be 
calculated by locating the position of the nearest singularity t c in the t plane. This singularity 
occurs when the modulus (3.147) becomes — oo, i.e. when 



e F (t c ) 



1 



(3.152) 



This is also the branch point in (3.146). It can be easily checked that this happens when 

tc = —■ 

3vr 



(3.153) 
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Figure 23: The Fermi level in the quartic potential with negative coupling t < 0. The nearest 
singularity corresponds to the critical value in which ep reaches the maximum of the potential. 



This has a nice interpretation in terms of the fermion picture. Since t c is negative, it corresponds 
to an inverted quartic potential and to the precise value of the parameter at which the Fermi 
level reaches the maximum of the potential, see Fig. 23. 

We can now try to calculate large N instanton effects in MQM. As in the case of standard QM 
with a quartic potential, we consider negative values of the coupling t = — k < 0. For this range 
of values, the vacuum at the origin is unstable and we should expect an instanton configuration 
mediating vacuum decay. In principle, one should write down an instanton solution with "small" 
action and calculate the path integral around it. This solution can be found by tunneling one 
single eigenvalue of the matrix M, which has an action of order 0(1), i.e. we consider the matrix 
instanton, 



where q c {t) is the bounce (3.21) with coupling A = g s . In principle, one could expand the 
path integral of MQM around this configuration and compute quantum planar fluctuations to 
determine the large N instanton action (this calculation was originally proposed in [104]). 

However, the fermion picture, which gives us a compact way of computing the planar ground 
state energy, should also give us an efficient way to compute the large N instanton action in a 
single strike. In this picture, the ground state is given by a filled Fermi level. As in any Fermi 
system, tunneling effects will first affect fermions which are near the Fermi surface. An instanton 
configuration with small action (i.e. of order O(N )) can then be obtained by tunneling a single 
fermion out of the N particles in the Fermi gas. Since at large N we can use semiclassical 
methods, the instanton action of such a fermion is just given by the standard WKB action, 



where a, b are the turning points associated to the Fermi energy e_p, and they are non-trivial 
functions of the 't Hooft parameter k. They are defined by the equations (3.146) with t = —k 
and an extra minus sign for b 2 . There is an extra factor of 2 due to the symmetry of the problem, 
and the factor of N is due to the fact that the effective Planck constant in this problem is 1/iV, as 
we remarked in (3.135) (we are setting H = 1). The integral in (3.155) can be explicitly computed 
by using elliptic functions, and the final result is 



M c (t) = diag (0, • • • ,0,g c (t), (),••■ 



0) 



(3.154) 




(3.155) 




(3.156) 
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where the elliptic modulus is now given by 



k 2 = (3.157) 

The above function has the following expansion around k = 0, 

, 4 , /16e\ 17k 2 125k 3 
AW=3-«log(-) + - nr + w + .... (3.158) 

This is precisely the expected structure (3.107) for a large N instanton action: the leading term 
is the action for the instanton (3.26) in the N = 1 quantum mechanical problem, the log term 
is a onedoop factor in disguise, and the rest of the series is a sum of loop corrections in the 
background of the "classical" instanton. An interesting property of A(n) is that it vanishes at 
the critical value 

/2 

K = t c = 2^-, (3.159) 

3"7T 

which is indeed the singularity (3.153) in the complex t plane signaling the convergence radius 
of the genus g ground state energies. It corresponds to the critical point at which the Fermi sea 
reaches the local maximum; at this point the action for tunneling must indeed vanish since the 
cndpoints in (3.155) collide: a(t c ) = b(t c ). 

Although we have just computed here the leading order term in the 1/N expansion of the 
ground state energy, there is a full series of higher genus corrections to this result S g (t), g > 1, 
which play the role of the genus corrections to the free energy discussed above. One would expect 
that these corrections display the asymptotic behavior (3.102), where A(t) is the action of the 
large N instanton calculated in (3.156). This was indeed verified in [98]. 

4. Non-perturbative effects in matrix models 
4.1 Matrix models at large N: General aspects 

In this subsection we review some basic facts about matrix models. For more detailed expositions 
and explanations, the reader is encouraged to look at [49, 92, 58]. 

We will consider matrix models for an N x N Hermitian matrix M, with a potential V(M). 
For the moment being we will assume that this is a polynomial potential, 

where the g p are coupling constants of the model. The partition function is defined by 

= »ra / W,M) . («) 

where g s is a coupling constant. Matrix models have a U(N) "gauge" symmetry 

M -»• UMU\ (4.3) 

therefore one can go to the "diagonal gauge" and write this partition function in terms of the 
eigenvalues of M, denoted by x« (this is very similar to the gauge (3.126) in matrix quantum 
mechanics). The resulting N- dimensional integral is given by 

i r N J 

Zl (N,g s ) = -LJ n^A 2 (*)e-^^> (4.4) 

T i=l 
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where A(x) is the Vandermonde determinant introduced in (3.131). Here 7 is a contour in the 
complex plane, which we take to be the same for the N eigenvalues, and which makes the integral 
convergent. Let us denote the critical points of V(x) by 

xi,- ■■ ,x d . (4.5) 

It is possible [59] to choose d integration contours 7^ in the complex plane, k = 1, ■ • ■ , d, going 
to infinity in directions where exp(— V{x)/g s ) decays exponentially, and in such a way that each 
of them passes through exactly one of the d critical points and is a steepest descent contour. 
In general, the contours depend on the argument of g s , but we can always make a change of 
variables so as to reabsorb this phase in the coefficients of V(x). The original contour 7 can be 
written as a linear combination of such contours, 



7 = E C ^- ( 4 - 6 ) 

fc=i 



We can then write 



Z 1 (N,g s )= C^---C^Z(N U ... ,N d ). (4.7) 

Nl +-+N d =N 



N 



(4.8) 



In this formula, 

Z(N 1 ,...,N d )= — 1 — / ... f f[^A 2 (A)e-^- y ( A * 

and we have split the eigenvalues in d sets, 

{^}i k =i,-,N k , k = l,--.,d. (4.9) 

Each integral (4.8) defines a possible "background" or vacuum of the original matrix model. 
There are three important remarks to be made about the above expansion: 

1. The backgrounds (4.8) of the matrix model are in one-to-one correspondence with gauge 
symmetry breaking patterns, 

U(N) -»• U(Ni) x • • • x U(N d ). (4.10) 

This is similar to the vacuum structure (3.58) of CS theory on lens spaces L(d, 1). We will 
denote these vacua by 

(N lr --,N d ). (4.11) 

2. These sectors appear as an artifact of the saddle-point approximation. The original in- 
tegral (4.4) is perfectly well-defined, and it is only our will to treat the integral in this 
approximation which led to the appearance of these sectors. 

3. The expression (4.7) is the analogue in random matrix theory of a trans-series solution. 
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A choice of background in the matrix model will define a choice of a perturbative sector. As 
we will see, the remaining sectors can then be regarded as instanton sectors w.r.t. the chosen 
background. A background of the form 

(TV, 0, ■■•,()) (4.12) 

is called a one- cut background, for reasons which will be clear in a moment, while the generic 
background will be called a multi-cut background. 

Let us now choose a fixed background and let us consider the integral (4.8) on this back- 
ground. It can be computed in a saddle-point expansion at small g s but keeping fixed the so-called 
partial 't Hooft couplings 

U = gsNt. (4.13) 

This means that we are doing a large N expansion of the matrix integral. The total 't Hooft 
parameter is 

d 

t = g s N = Y,Ni. (4.14) 

i=i 

One can use the standard large N counting arguments to see that the structure of the free energy 
is of the form 

oo 

logZ(M,-" ,N d ) = Y J 9 2 s g - 2 F g {t l ,--- ,t d ). (4.15) 

9=0 

In the one-cut case d = 1 this is the standard 1/N counting briefly reviewed in the previous 
section. In the multi-cut case one has to do a slightly refined analysis, see the Appendix in the 
published version of [27] for an explicit derivation. 

Example 4.1. The Gaussian matrix model. The Gaussian matrix model is defined by the matrix 
integral (4.4) with the potential 

V(x) = ^. (4.16) 

If g s > 0, the integration contour is simply 7 = R, the real axis. We will denote the partition 
function of this model by Z G (N, g s ). This is one of the few cases in which the matrix integral 
can be computed exactly at finite N, 

N 2 /2 

Z G (N,g s ) = ^L_G 2 (tf + l), (4.17) 
where G 2 (N + 1) is the Barnes function 

AT-1 

G 2 (N + 1) = JJ i\. (4.18) 

i=0 

In this case, the large N expansion (4.15) follows from the asymptotics of this function. One 
finds, 

i^) = ^ 2 (logt-^), 

*f (*) = ~ log t + ^ log g s + C'(-l), (4.19) 
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where -E>2# are Bernoulli numbers. Notice that (t) depends also on g s , but usually this piece 
(as well as the constant involving the zeta function) is not taken into account. 

The genus g free energies F g (ti, ■ ■ ■ ,t d ) appearing in (4.15) can be expanded around t\ = 0. 
The first terms in this expansion are just the Gaussian free energies for the different Vs. Once 
these terms are subtracted, the resulting quantities 

d 

F g (t lr -- , td )-Y,Ff(t t ) (4.20) 
i=i 

are analytic at the origin, in agreement with the general arguments put forward in the previous 
section. We will regard (4.15) as our perturbative expansion. It is interesting to note that 
this expansion is a generalization/deformation of the standard saddle-point expansion of one- 
dimensional integrals. Indeed, let us write (4.8) as 

We want to take a limit where g s is small and t is fixed, so that ./V is large. The two terms inside 
the parenthesis are then of the same order (i.e. O(N)) and it is clear that the 't Hooft parameter 
controls the strength of the Vandermonde interaction. Let us suppose that t is very small, so 
that we can neglect this interaction. In this limit (4.8) factorizes into a product of standard 
saddle-point integrals, 

d 

Z(N 1} --- ,N d )^H(f lt (g s )) N \ t^O, (4.22) 
i=i 

where 

fM = l^e~^*e-^ (4.23) 

and we approximated this integral by a saddle-point expansion around X{. When t is no longer 
small, we have to take into account the Vandermonde determinant. Since this induces a repulsion 
between eigenvalues, they will no longer sit at the saddle points of V(x): the Nk eigenvalues in 
the k-th set will sit at an interval or arc Ck around the k-th extremum. As we will see in a 
moment, when N is large but the 't Hooft parameters tj are fixed, these arcs are compact, 
and the distribution of the eigenvalues on these arcs is given by a density function p(A), whose 
support is the union of the intervals. The problem of finding this equilibrium distribution can 
then be regarded as a deformation of the saddle-point technique for standard integrals, where 
the deformation parameters are the 't Hooft parameters tj. 

4.2 The one-cut solution 

The determination of the free energies F g (ti) and of the density function p(X), for a given matrix 
model potential and background, has a long story, which starts in [28] and culminates in [57]. We 
will now review some of the relevant results, which we will need to develop the instanton calculus 
in matrix models. To begin with, we will consider as our background a one-cut background. 
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In this case, there is a single 't Hooft parameter and the free energy has a perturbative genus 
expansion of the form 

oo 

F = Y J F g {t)9 2 r 2 . (4-24) 

3=0 

Another important set of quantities in a matrix model are the connected correlation functions 

W h (p 1 ,..., Ph ) = (Tr— i — ...Tr—!— ) , (4.25) 
\ pi — M p h - M I (c) 

where the subscript (c) means connected. These correlation functions are generating functions 
for multi-trace correlators of the form 

W h ( Pl , ..., Ph )=J2 Wl+ i 1 n h+ i ^ M ' ni ■ ■ ■ Tr M«-) (c) , (4.26) 
m>iPi '"Ph 

and they further have a g s expansion of the form 

oo 

W fc (pi, . . . ,p h ) = Y,9s g+h ~ 2 W 9th ( Pl , . . . ,p h ). (4.27) 

9=0 

At large N the one-cut background is characterized by a density of eigenvalues p(X) which 
has support on a single, connected interval C = [a, b] in the complex plane. This density is 
completely determined by the condition that the so-called effective potential on an eigenvalue, 

V cS (X) = V(X)-2t J dAXAOloglA-A'l, (4.28) 

has to be constant -at fixed 't Hooft coupling- on the interval C: 

V e g(\) = m), AGC. (4.29) 

A quantity which is closely related to the density of eigenvalues is the planar resolvent, which 
is nothing but the quantity Wo t \(p)/t, where Wq,i(p) has been introduced in (4.27). The planar 
resolvent can be computed as 

UH>(p) = J dA^. (4.30) 
It also satisfies the asymptotic condition 

uj (p) ~ -, p -)• oo, (4.31) 
p 

which follows from the normalization of the density function 

d\p(\) = l. (4.32) 



IC 

Once the resolvent is known, the eigenvalue density follows as 



p(A) = --L(wo(A + ie) -wo(A-ie)). (4.33) 

Z7T1 
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It can be seen that the condition (4.29) determines the resolvent as 

1 / dzV'(z) f(p-a)(p-b)\* , A 

Here the integration is around a closed contour which encircles the cut C. By deforming the 
integration contour, this solution can also be written as, 

"o(p) = Y t (V(p)-y( P )), (4.35) 

where y(p) is a function on the complex plane which has a branch cut along C, called the spectral 
curve of the matrix model. In the one-cut case, the spectral curve has the structure 

y{p) = M(p)y/(p-a)(p-b), (4.36) 

where M(p), known as the moment function, is given by 

„(„, -f*Lm , 1 (4.37) 

with the contour of integration being around the point at oo. The endpoints of the cut follow 
from the asymptotic behavior of the resolvent (4.31), leading to the equations 



i 



Az v ' {z) - 0, 



Jc 27ri ^( z - a )( z -b) 

It turns out, with the exception of the free energies at genus zero and one, and the planar 
one-point function Wo,i(p), the quantities F g (t) and W 9j h(pi, ■ ■ ■ ,Ph) c & n be computed in terms 
of the spectral curve alone. More precisely, knowledge of the endpoints of the cut, a and b, and 
of the moment function, is all one needs in order to compute them. This was first made clear 
in [28, 9, 8] and later culminated in the geometric formalism of [54, 36, 57]. For example, the 
two-point correlator at genus zero is given by [9] 

Wo2ip , „ _ _JL_ ( «-}(■> + «)(. + ») + <* _ A (4 39) 

and only depends on the end-points of the cut. On the other hand, the genus-zero free energy, 
Fo(t) is given by 

F (t) = -\! d\p{\)V(\) - \t 2 m- (4-40) 

The derivatives of the planar free energy can be computed in terms of the effective potential. 
One finds (see for example [77]) 



d t F (t) = -tm = -v cS (b), 

d?F (t) = -d t V cS (b) = 2log 



b- a Ml) 
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4.3 The multi-cut solution 

We will now consider the more general (and more difficult) background: the multi-cut solution. 
In this case the support of the eigenvalue distribution is a disjoint union of n intervals 

n 

C = {JC t . (4.42) 
i=\ 

The density p(A) now satisfies the equation, 

V eS (X) =Ti, AG d, (4.43) 

where Tj are constants in each interval Cj. This generalizes (4.29) to the multi-cut case. 

The way to implement the multi-cut solution at the level of the planar resolvent is to require 
wq(p) to have 2n branch points. The solution now reads, 



This can be also written as in (4.35), in terms of a spectral curve, which is now hyperelliptic, 



2n 



y(p) = M (p) Yl Vp^k. (4.45) 



k=i 



Here, M{jp) is the multi-cut moment function, which is given by the obvious generalization of 
(4.37) to the multi-cut case. In order to satisfy the asymptotics (4.31) the following conditions 
must hold: 

z e V'(z) 



Sen = ^f ^=7 - £ = 0,l,---,n. (4.46) 

In contrast to the one-cut case, these are only n + 1 conditions for the 2n variables xu representing 
the endpoints of the cut. The remaining n — 1 conditions are obtained by fixing the values of 
partial 't Hooft parameters: 

t i = g s ^d\p(X). (4.47) 

The planar free energy satisfies in addition the equation, 

dF 8F 



dU dt i+ i 



= r m - r i5 (4.48) 



where are the quantities appearing in (4.43). 

We can write the multi-cut solution in a very elegant way by using contour integrals of the 
spectral curve. First, the partial 't Hooft parameters are given by 

J d 

We now introduce dual cycles T>i cycles, i = 1, ■ ■ ■ ,n — 1, going from the Cj + i cycle to the Ci 
cycle counterclockwise, see Fig. 24. In terms of these, we can write (4.48) as 

dF 8F If 
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Figure 24: A two-cut spectral curve, showing two contours C1.2 around the cuts where eigenvalues 
sit. The "dual" cycle T>\ goes from C 2 to C\. 

4.4 Large N instantons and eigenvalue tunneling 

In the sum (4.8) over possible "vacua" of the matrix model, any two vacua are related by a 
redistribution of the eigenvalues, i.e. by a sequence of operations of the form 

Ni -> Ni - 1, Nj -»• Nj + 1. (4.51) 

This can be interpreted as a process in which an eigenvalue leaves the cut around the i-th critical 
point and "tunnels" to the cut around the j-th critical point. Therefore, given a choice of 
background or perturbative vacuum, all the other vacua in (4.8) (which are then regarded as 
instanton sectors) can be obtained from the reference background by eigenvalue tunneling. The 
fact that instantons in the matrix model are due to eigenvalue tunneling was first pointed out in 
[44, 112]. 

Notice that, when an eigenvalue tunnels as in (4.51), the r.h.s. of (4.22) changes, at leading 
order, by 

exp f-l(y( Xj ) - V( Xi ))]. (4.52) 
We can think about this as the action of the classical instanton connecting the two vacua: 



(••• ,N h --- ,AT jV ..)^(... , JV< — 1, - - ■ ,Nj + l,---). (4.53) 

As we have explained in the previous section, there is a large N instanton built upon this classical 
configuration. The classical action (4.52) gets corrected to 

exp (-^-(VesiXj) - VesiXi))) (4.54) 

where the effective potential is defined in (4.28). This can be regarded as the action of the large 
N instanton connecting the vacua. It follows from (4.48) and (4.43) that the large N instanton 
action can be computed in terms of the genus zero free energy as 

V eS ( Xj ) - V cS ( Xl ) = -^--^. (4.55) 

When t — >• 0, we recover the classical instanton action (4.52). The picture of eigenvalue tunneling 
is still valid for finite t, but now the tunneling is between two cuts rather than between two 
saddle points. A graphical depiction of an ^-eigenvalue tunneling in a cubic potential, from the 
background (N, 0) to the background (N — £,£), is shown in Fig. 25. 
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Figure 25: Eigenvalue tunneling in a cubic matrix model. 
4.5 Large N instantons in the one-cut matrix model 

A particularly simple case of eigenvalue tunneling occurs when the reference background is the 
one-cut configuration. The elementary instanton configurations correspond in this case to tun- 
nelings of the form 

(N, 0, • • • , 0) — >• (JV — 1, 0, • • • , 0, 1, 0, • • • , 0). (4.56) 

The large ./V instanton action depends then only on t = g s N and it has the structure discussed 
in the previous section. Since this is the simplest case, we will now explain how to calculate 
instanton contributions when the reference background is a one-cut solution. It turns out that 
this computation can be done in at least three different ways. The first method is based on 
a direct calculation of the matrix integral. It was introduced by F. David in [44] and further 
clarified in [75, 77]. The calculation of [44, 75, 77] is done for the double-scaled matrix model, in 
which the 't Hooft parameter is near its critical value. The correct calculation of the instanton 
contribution at generic values of the 't Hooft parameter was presented in [100]. There is a 
second method, which regards the instanton sectors in the one-cut matrix model as limits of the 
generic multi-cut vacuum. This method has the advantage of giving general expressions for the 
k-th. instanton sector, and it was presented in [101]. Finally, a third method was introduced in 
[94], based on the method of orthogonal polynomials, which makes it possible to calculate the 
instanton sectors as trans-series solutions to difference equations, and therefore it mimicks the 
structure developed for ODEs. We will present the three methods in turn, skipping some of the 
details which can be found in the original references. 



4.5.1 A direct calculation 

We thus consider a one-cut matrix model in which the effective potential has the form depicted in 
Fig. 26: it is constant along the location of the cut C = [a, b], and there is a critical point xq which 
corresponds to another possible vacuum. We will consider the instanton sectors corresponding 
to eigenvalue tunneling from the cut at C to the critical point xq. These instanton sectors are 
labelled by a positive integer £, which is the number of tunneling eigenvalues. The total partition 
function, summing over all these instanton sectors, is given by 

Z(N) = Z (0 \N) ^1 + C e ZW(N)j , (4.57) 

where for convenience we have implicitly normalized the sectors with £ ^ by the perturbative, 
one-cut partition function Z^ \N), and C is a constant. Notice that this total partition function 
has the trans-series structure considered in section 2. The free energy also has a trans-series 
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M(x ) = 




Figure 26: The effective potential for the one-cut matrix model. The point xq is a critical point where 
eigenvalues can tunnel to, leading to instanton configurations. 



structure, 

/ oo \ oo 

F(N) = F<-°\N) + log ll + ^2c e zW(N) \ =Y,C e F^(N). (4.58) 



The instanton partition functions (N) can be computed directly in terms of the defining 
matrix integral: we consider an integral in which N — £ eigenvalues are located in the cut C, 
while for the t eigenvalues which have tunneled, the integration is made around a steepest- 
descent contour passing through the critical point. In particular, for the one-instanton sector 
one finds [75] 

,r . f . N-l 

Z<X>(N) = / d*e-£ v « / TTdA,A 2 (x,A 1 ,...,A iV _ 1 )e-^^ 1 ^ ) , 

(4.59) 

where the first integral in x is around the critical point at xq, along a saddle-point contour which 
we have denoted by X. The rest of the N — 1 eigenvalues are integrated around the contour Xq 
corresponding to the background one-cut configuration. The overall factor of N in front of the 
integral is a symmetry factor, counting the N possible distinct ways of choosing one eigenvalue 
out of a set of N. One can easily write similar integrals for the ^-instanton contribution, but 
these will be easier to calculate with the second method below. The integral (4.59) can be written 
as 

z<1,(w »=^^vrL d ^ w ' (4 - 60 » 

where 

f{x) = <det(xl - M') 2 }^ e-^ V{x) . (4.61) 

The notation in these equations is as follows: Z^°\N) is the "perturbative" partition function 
in the one-cut background with iV eigenvalues, the average is defined as 

\ U /N — ~ M i ^JV vn \ V 1 - *) 
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and it is calculated again in the one-cut bacgrkound. Finally, M' is an (TV— 1) x {N— 1) hermitian 
matrix. We conclude that the one-instanton sector can be evaluated by calculating correlators 
in the perturbative sector. In fact, by making use of the familiar relation 



det(xl — M) = exp (trln(xl - M)) 



we obtain 



(det(xl - M) 2 ) = exp 



X^((trln(*l-M)r> 



(c) 



s=l 



(4.63) 



(4.64) 



which is written in terms of connected correlation functions. The correlation functions appear- 
ing in (4.64) are nothing but integrated versions of the Wh correlators in (4.25), evaluated at 
coincident points. Let us define 



A, 



dpi-- - / dp h Wg th (p 1 ,--- ,p h ) 



x\=—=x h =x 



La J 2 n ~ 2k + l 
A n (x;t) = X] ( n _ 2fc + 1)! A k,n-2k+l{x;t), H > 1. 



(4.65) 



k=0 



In this notation, the general perturbative formula for the determinant reads 

<det(xl - M) 2 ) = exp (f^g^M^t)) . 



\n=0 



where A n (x;t) is the n-loop contribution. We have, for example, 

A {x;t) = 2A 0:1 (x;t), Ai(x;t) = 2A , 2 (.x;t). 



(4.66) 



(4.67) 



The integration constants involved in the integrations in (4.65) may be simply fixed by looking 
at the large x expansion of the correlators. Next, we define the holomorphic effective potential, 
which combines the matrix model potential together with Ao(x;t), as 

H,eff(x; t) = V{x) -2tJ dpuoo(p) = V(x) -2tJ dp p(p) log(x - p). (4.68) 



It satisfies 
as well as 



ReV hjeS (x;t) = V eS (x), 
where V e g(x) was earlier defined in (4.28). Altogether, one finally has for the integrand 



f(x) = exp (--F hiCff (x;0 + f>r^4n(a;;0 J , 
V 9s n=l / 



where 



t' = g s (N-l) =t-g s . 



(4.69) 
(4.70) 

(4.71) 
(4.72) 



This shift in the 't Hooft parameter is due to the fact that the correlation function involved in 
(4.61) is computed in a matrix model with N — 1 eigenvalues (recall we removed one eigenvalue 
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from the cut). Since we are computing the one-instanton contribution in the theory with N 
eigenvalues, we thus have to expand (4.71) around t. This gives further corrections in g s , and 
one finds 

f(x) = exp (-jV^(x) + $(x)\ , (4.73) 



with 

oo oo 

= EsT 1 *»(*) = £ sT 1 

n=l 

One has, for example, 



n=l 



("I) 



n-1 



n-1 



rv. 



fc=0 



$i(x) =Ai(x) + dtV h ,eB(x). 



(4.74) 



(4.75) 



Further explicit results for these quantities can be found in [100]. In the expression (4.73) all 
quantities now depend on the standard 't Hooft parameter t for the model with N eigenvalues, 
and we have thus dropped the explicit dependence on t. The derivatives with respect to t can be 
performed by using standard results in one-cut matrix models, some of which were listed above. 
One may now proceed with the integration of f(x), 



dx exp 



1 



V h , eS {x) + <S>(x) 



(4.76) 



If we wish to evaluate this integral as a perturbative expansion when g s is small, we can do it 
using a saddle-point evaluation [75, 77]. The saddle-point condition 



(4.77) 



requires xo to be a critical point of the effective potential, as anticipated above. If we use the 
explicit form of the spectral curve (4.36) we find the equivalent condition 




M(x ) = 



Figure 27: The spectral curve y(x) has a singular point at the nontrivial saddle xq. 



M(x ) = 0. 



(4.78) 



Geometrically, the spectral curve is a curve of genus zero pinched at xq, as shown in Fig. 27. 
This was observed in [111] in the context of spectral curves for double-scaled matrix models, and 
their relation with minimal strings (see also [83]). Of course, it can happen that there are many 
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different critical points of the effective potential, therefore more than one solution to (4.78). In 
this case, there will be various instantons and we will have to add up their contributions (the 
leading contribution arising from the instanton with the smallest action, in absolute value). 

The calculation of (4.76) is now completely standard, and it reduces to Gaussian integrations. 
The result is 

L di/w =\Z^ exp (-^ h ' tfW+#i(xo) ) ( l+ § 9 "4 (4 - 79) 

where the f n can be systematically computed in terms of the functions $> n (x) and their deriva- 
tives, evaluated at the saddle-point xq. An explicit expression for /2 can be found in [100]. 

With these ingredients, we can already calculate the one-instanton contribution to the free 
energy. Using (4.58), we find 

^^^"'w^wL^ 14 (4 - 80) 

The quotient of perturbative partition functions is easy to calculate in terms of the free energies, 



Z(°)(N) 

where 



Z( °X m l) = exP (m ~ F(t)) = exp (f^n) , 



(4.81) 



\n=0 



One has, for example, 



L 2 J (_-\\n-2k+l 

e^E ;j 2k + lv ^ M mi (4. 8 2) 

Go = -d t F (t), Qx = \ d?F (t). (4.83) 
Putting together (4.79) and (4.82) above, we finally find that has the structure 



where, up to one loop, we have: 

A = V hjeS (x ) - Go{t), 



(4.85) 



We can give explicit expressions for all the quantities involved in (4.84) in terms of data 
which depend only on the spectral curve (4.36). First of all, by using (4.83), (4.41) and (4.69) 
we find 

A = V hie ${x )-V h , eS (b) = ^ dzy(z), (4.86) 

which is the instanton action (here, we use the fact that Vh e ff(&) = V e s(b)). As pointed out in 
[111], this expression also has a geometric interpretation as the contour integral of the one-form 
y(z) dz, from the endpoint of the cut C to the singular point xq (see Fig. 27 and Fig. 28, left 
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side). To compute up to one-loop, we must compute $i(x), given in (4.75). One can find 
the result for Ao t 2(x;t) (which enters in the expression of Ai) simply by integrating the first 
formula in (4.39) '[77] 



A , 2 {x;t) = log 1 + 



Using that (see for example [49]) 



one further finds, 



and one then obtains 



d(tuj (p)) 
dt 



dtV hte fi(x) = -4 log 



x-(a + b)/2 
y/(x -a)(x- b) 

1 

V \P - a){p - b)' 
\/x — a + \J x — b 



-log 2. 



4 log 2, 



&i(x) = — log (x — a)(x — b) 



(4.87) 



(4.88) 



(4.89) 



(4.90) 



Adding to <3?i(x) the result for Qi{t), which follows from (4.41), it is simple to put all expressions 
together and obtain the contribution, F^i, of the one-loop fluctuations around the one-instanton 
configuration, 

(4.91) 



F u = -i 



\ 2ttM'(x ) (x - a)(x - b) 



_, 5 • 

2 



This formula is valid for any one-cut matrix model where the potential has an extra critical point 
xo- Notice that if xo is a local maximum of V e s(x), one will have that M'(xo) < 0, and hence 
Fi t i will be real. This result extends previous calculations in [44, 75] to one-cut matrix models 
away from the critical point. 



Example 4.2. Large N instanton in the quartic matrix model. The quartic matrix model is 
defined by the potential 

V(x) = \x 2 - ^x\ (4.92) 

where we follow the normalization of [94]. This potential has three critical values, namely 

2 



x = 0, 



x = ± 



(4.93) 



so the generic matrix model based on this potential has three cuts. We can however consider 
the one-cut model where the eigenvalues sit around the origin x = 0. The spectral curve for this 
one-cut background can be easily obtained with the standard matrix model techniques reviewed 
above. It is given by 

y = M{x)\/x 2 - 4a 2 

where the moment function is 



M(x) = 1 



^2 ^2 

-a x 



3 6 

The cut is located at C = [—2a, 2a], and the endpoints are determined by 



(4.94) 
(4.95) 



a 



X 



i - Vi-tx 



(4.96) 
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The moment function has two zeros which give two non-trivial saddle-points of the effective 
potential, namely ±xo with 

Ul- X a 2 ) . (4.97) 



A 

As t — > 0, they become the other critical points of the potential, 

2 

x -> ±-=. (4.98) 
v A 

Since the potential is symmetric, there are too one-instanton solutions, corresponding to one 
eigenvalue tunneling from C to the two saddles ±xo- Both instantons have the same action, 
which is computed by integrating the spectral curve as in (4.86), 



A V 2 / \ 6 



- ^ ( a *A - 4) (log (£) - 21og ^1 - 2* + - ^ 
This has the small t expansion 

l . / . . /AAA 



(4.99) 



A{t) = 3X +t {~ l + lQg V 12 J ) + ° {t ) - (4 - 100) 
The first term (i.e. the instanton action as t — > 0) is given by 

V (-j=) ~ V(0), (4.101) 

in agreement with (4.52). The structure of the next-todeading correction is in agreement with 
the general expectation (3.107). 

Example 4.3. Large N instanton in the cubic matrix model. Let us consider the matrix model 
analogue of the Airy integral (2.62). The potential is 

x 3 



V(x) = -e 1K x + — . (4.102) 



This has two critical points 

where we have introduced the variable 



x l,r = T £i/2 ( 4 103 ) 



C = e iK . (4.104) 



The most general background of a matrix model based on this potential is a two-cut configuration, 
labelled by (N±, N2), where N±, N2 are the number of eigenvalues near the critical points x L,R , 
respectively. Let us consider the one-cut background where all the eigenvalues sit near x L , i.e. 

(N 1 ,N 2 ) = (N,0). (4.105) 

The endpoints of the cut (a, 6) are determined by the equations 

x (4-0=t, 5 2 = 2(C-x 2 ), (4.106) 
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where 

a = - XQ + S, b = -x -S. (4.107) 

The spectral curve is given by 

y [x) = {x- x )^x 2 + 2xx + 3x1 - 2 C- (4.108) 
The effective potential is (up to a constant) the integral of y(x), which gives 

V cS (x) = \ (x(x - s ) " 2C) J x 2 + 2xx + 3x 2 - 2( 

-l„,4-o, g (, + , + ^^; 

The instanton action is given by 

A = J y{x)dx = V cfi {xo)-V c{i {b), (4.110) 
and it has the small t expansion 



(4.109) 



A = -^C 3/2 + t (-1 + log ^ ) + 0(t 2 ), (4.111) 



and as expected the leading term as t — > is 

4 



3 C 3/2 = V(x R )-V(x L ). (4.112) 
4.5.2 Large iV instantons from multi-cuts 

It is clear that the instanton configurations in the one-cut matrix model that we have just 
analyzed, in which I eigenvalues tunnel to other critical point, can be regarded as particular 
cases of a two-cut configuration with 

N 1= N-£, N 2 = £, I -v. A. (4.113) 

(£) 

Therefore, we should be able to calculate the general instanton partition function Z N as a 
particular case of a two-cut partition function, i.e. we expect 

ZW(N) = Z(N -£,£). (4.114) 

This was made explicit in [101] and gives compact formulae for the general ^-instanton sector of 
the one-cut matrix model. We can evaluate the r.h.s. of (4.114) by expanding the free energy of 
the two-cut matrix model 

F (ti = t - g s £,t 2 = g s £) (4.115) 

in powers of g s . There is however one subtlety: the free energies F g (ti,t2) are not analytic at 
t 2 = 0. Geometrically, this corresponds to the fact that we are expanding around a configuration 
in which the second cut of the spectral curve is completely pinched, as shown in Fig. 27. However, 
as we mentioned in (4.20), this non-analyticity is due only to the Gaussian contribution to the 
two-cut integral. Therefore, the functions F g (ti,t2) defined by 

F g (h,t 2 ) = Ff(t 2 ) + F g (ti,t 2 ), (4.116) 
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where F g (t) are the genus g Gaussian free energies written down in (4.19), are analytic at £2 = 0. 

Physically, the reason for the appearance of this singularity is that in this problem t2 = ig s , 
and I is small as compared to N. Therefore, it is not appropriate to treat the integration over 
the t tunneling eigenvalues from the point of view of the large N expansion. Instead, they should 
be integrated exactly. This argument also suggests that, in order to regularize the computation, 
we should subtract F G (t2) from the total free energy and at the same time multiply by 
the exact partition function Zf, which is given in (4.17). The appropriate expression for the 
partition function around the ^-instanton configuration is then 



(F 9 (t-e 9s ,£g s )-F g (t)) 



= Zf exp 

L <y>o 

By expanding in g s , this expression leads to the general formula 



(4.117) 



Y\(2 ft+mi -2) (4.118) 

x 2^ 2^ 2^ k\ mi \ . . . m k \ 91 " 9k { ' 

In this equation we have introduced the following notations: the m-th derivative of F g is taken 
w.r.t. the variable s, which is defined as 

s = \{t 1 -t 2 ). (4.119) 

All derivatives are evaluated at t\ = t and t2 = 0. The exponential factor involves the instanton 
action, which is given by 

A(t) = d s F (4.120) 

and 

q = exp(d 2 s F ). (4.121) 

At leading order in g s , we have 

E,(2< ?l +m I -2) 

2^2^ 2^ k\m l \...m k \ 91 " 9k { ' 

k mi>0 gi>1 ^ ' (4.122) 

= 1 - gsildsF^t) + j d 3 s F (t)) + 0(g 2 s ). 

By using the general formulae for the free energies in the two-cut model, one can check [101] that 
(4.118) reproduces the one- instanton contribution to the free energy computed in the previous 
subsection and generalizes it to arbitrary t. 

4.5.3 Large iV instantons from orthogonal polynomials 

So far we have analyzed the "perturbative" sector of the one-cut matrix model from the point of 
view of the spectral curve: the procedure explained in [36, 57] gives a systematic way to compute 
the "perturbative" 1/N corrections in terms of data of this curve, and the results reviewed 
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above show that the formal trans-series giving the instanton corrections can be also computed 
using data of the spectral curve. For example, the instanton action is given by a period of the 
differential y(x)dx along a cycle, and the loop corrections can be computed from the knowledge 
of the curve only. 

However, in the one-cut case there is a completely different way to obtain the 1/N cor- 
rections in the perturbative sector, by using orthogonal polynomials [24, 25] (see [6] for some 
recent developments). It was shown in [94] that this procedure can be generalized to compute 
systematically multi-instanton effects in the one-cut sector (partial results in this direction can 
be found in [67, 4, 113]). This technique is very powerful in order to go to higher loops in sim- 
ple matrix models. In addition, it is formally very similar to the trans-series method that we 
developed in section 2 and therefore it makes possible to apply many of the techniques typical 
of the theory of resurgence (the connection to resurgence theory was further developed in [10]). 
Another advantage of this method is that it can be applied to unitary matrix models. 

The basic idea of the technique of orthogonal polynomials is the following (see [49, 92] for 
reviews): in the matrix integral (4.4), if we regard 

d ^ = e -^(A)d^ (4123) 

as a measure in R, we can introduce orthogonal polynomials p n (X) defined by 

d/ip n (\)p m (\) = h n 5 nm , n > 0, (4.124) 



where p n (X) are normalized by requiring the behavior p n (X) = X n + • • • . One then easily finds, 

N-l N 



Z= Y[h i = h$l[rf- i , (4.125) 

i=0 i=l 

where we have introduced the coefficients 

r k = k>l, (4.126) 

which appear in the recursion relations for the p n (X), 

(A + s n )p n (X) = Pn +i(X) + r ra p„_i(A). (4.127) 

In this subsection F will denote the normalized free energy, which is obtained by subtracting the 
Gaussian free energy, 

F = logZ-logZ G . (4.128) 
At finite N, F is given by the following formula: 

u k=l 

where is the coefficient ho for the Gaussian model. In order to proceed, we introduce a 
continuous variable as iV — > oo, 

g s k^ z, 0<z<t, (4.130) 
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and we assume that in this continuum, N — > oo limit, becomes a function of z and g s , 

r k ^R(z,g s ). (4.131) 

It will be useful to consider the function 

R(z,9s) 



(4.132) 



which can be regarded as the continuum limit of rh/(kg s ). It is easy to see that, for polynomial 
potentials of the form 

V(M) = ^M 2 + ■■■ , (4.133) 

one has ~ kg s + - ■ ■ , therefore the function log(rk/(kg s )) is regular at k = and we can use the 
standard Euler-Maclaurin summation formula to evaluate (4.129). One then obtains [24, 25]: 



rt 00 

g 2 s F= / dz(t-z)\ogE(z)+J29. 
J o p=i 



2 P D iv 



(2p)\ dz^- 1 



(t - z) logZ(z,g s 



z=t 

2=0 



+ 



tg s 



(4.134) 



From this expression one can deduce the difference equation, 

F(t + g s ) + F(t- g s ) - 2F(t) = log H, 
which can be also obtained by starting from the identity 

Zn+iZn-i 



7 2 



r N . 



(4.135) 



(4.136) 



In order to compute the g s expansion of the free energy, one has to first find an expansion 
for the function R(z,g s ) of the form 



R^(z,g s ) = f29 2 s S Ro,2sU 



(4.137) 



s=0 



Once this expansion is plugged in E(z,g s ) and then in (4.134), the genus expansion follows. 
In order to obtain (4.137) one has to use the so-called pre-string equation. This is a difference 
equation for R(z, g s ) which can be derived as the continuum limit of the recursion relations obeyed 
by the coefficients (4.126). The pre-string equation can be explicitly written for any polynomial 
potential [25, 49]. For example, in the case of the quartic matrix model with potential (4.92), 
the difference equation for R(z,g s ) reads as 



R(z,9s){l ~ Y^{R(z,g s ) + R(z + g s ,g s ) + R(z - g s ,g s )} 



z. 



(4.138) 



These types of equations have a solution of the form (4.137), and they determine Ro jS (z) in terms 
of the Ro yS /(z), s' < s. When this solution is plugged in (4.135), one obtains the perturbative 
expansion of the total free energy in powers of g s , which is the standard 1/N expansion of the 
matrix model [24, 25]. 
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Difference equations, just like differential equations, admit trans-series solutions, and one 
could guess that the trans-series solution to the difference equation governing R(z,g s ) encodes 
the multi-instanton amplitudes of the full matrix model. This was first proposed in [94]. It was 
then verified that the results obtained with this technique agree with the previous techniques in 
the quartic matrix model [94, 10] and in the cubic matrix model [101] . To obtain the trans-series 
solutions to the pre-string equation, we consider a more general ansatz than (4.137), 



R(z,g s ) = Y,C e R {£) (z,g s ), (4.139) 
where R(°\z,g s ) is given by (4.137), and for £ > 1 we have 

oo 

R^(z,g s ) = e- eA ^R eA (z)(l + J>^, n+1 (z)), £ > 1. (4.140) 



n=l 



Once this ansatz is plugged in the difference equation for R(z,g s ), one obtains a recursive system 
of equations for the different quantities involved. The quantity A(z), which is a parameter- 
dependent instanton action, is determined by an equation of the form 

A'(z) = f(R , (z)), (4.141) 

where / is a function fixed by the difference equation. For £ = 1, n > 0, one obtains an equation 
which determines 

dRi >n (z) 



dz 



(4.142) 



in terms of Ri, n '(z) with n' < n. For n = 1, we have a differential equation for the logarithmic 
derivative, i.e. for 

1 dR hl (z) 



Ri,i{z) dz 



(4.143) 



The integration constant for R± t i(z) can be reabsorbed in the parameter C, and for A(z) and 
the Ri ;Tl (z), n > 1 the integration constants are fixed by using appropriate boundary conditions. 
For £ > 1, the difference equation determines in terms of Re, n ' with £ < £'. 

The trans-series structure of R(z,g s ) leads to a trans-series structure for the full free energy, 
as in (4.58). We will write it as 

oo 

F(t,g s ) = Y J C e F^(t,g s ), (4.144) 
e=o 

where 

oo 

F^(z,g s ) = e-WVO'Fviz)^ + ^ftn+iW). * > L ( 4 - 145 ) 

n=l 

Once (4.139) is known, one can plug it in (4.135) to deduce the F^\t, g s ), see [94] for more 
details. The resulting amplitude is nothing but the ^-instanton amplitude of the full matrix 
model. 
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Example 4.4. The quartic matrix model. As an example, let us present some results for multi- 
instanton corrections in the quartic matrix model with the potential (4.92). The perturbative 
solution (4.137) has been much studied since it was first worked out in the pioneering papers 
[24, 25]. The planar part is given by 

Ro,o(z) = |(l-VT^Ai). (4.146) 
As already noticed in [25], it turns out to be useful to express all results in terms of 

r = Ro,o(z). (4.147) 

For the higher g s corrections one finds, 



4 : 



3 (2-Ar) 

„ , \ 28A 4 r(5 + Ar) /A A . 

4A 6 r (111AV + 5728Ar + 7700) 
RofiW = ^ (2-Ar)" ' 

and so on. If we now plug the trans-series ansatz (4.139) in the difference equation (4.138), we 
find a system of recursive difference equations for the R^\z,g s ) which can be solved by using 
the ansatz (4.140). Let us focus on I = 1, the one-instanton solution. The first thing to compute 
is A(z), which corresponds physically to the instanton action. One finds, at leading order in g s , 

cosh(A'(z)) = 2^^. (4.149) 

This can be integrated to find A(z) up to an additive constant and an overall sign (since cosh z is 
even). Since z stands here for the 't Hooft parameter, both ambiguities can be fixed by requiring 
that, as z — > 0, the instanton action becomes (4.52). The result is 

Ar \ 



A(z) = - J dr cosh^ 1 (2^-^) (l 

= I r (Ar - 4) cosh" 1 ^--2j+ -L ^3(2 - Ar)(6 - Ar). 



(4.150) 



It can be checked that (4.150) coincides with the instanton action of the quartic matrix model 
computed in terms of its spectral curve in (4.99). Once the instanton action is known, we can 
proceed to compute R± ; i(z). The equation one obtains at the next order in g s is 

j^ = --coth(A>(z))A"(z), (4.151) 
which can be immediately integrated as 

/ \ -1/2 

R hl (z) = (smh(A'(z))) . (4.152) 

The rest of the coefficients can be found by integrating the resulting equations for Ri >n (z), 
and from (4.135) one finds the loop expansion of the one-instanton free energy, see [94, 10] for 
explicit formulae. The results agree with those obtained with the method of direct calculation 
in subsection 4.5.1 and with the result obtained from the multicut matrix model in subsection 
4.5.2. 
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4.6 Large N instantons, large order behavior and the spectral curve 

In the previous section we have discussed various techniques to compute large N instanton 
effects in one-cut matrix models. The total free energy, including multi-instanton sectors, is 
the analogue of the trans-series solution in this situation. We should then expect a connection 
between the large order behavior of the perturbative genus expansion and the instanton trans- 
series, generalizing the connections that we have seen in ODEs and in QM. Let us consider the 
total free energy 

F(g s )=g 2 F(g s ), (4.153) 

which is defined in such a way that the perturbative sector has no negative powers of g s . The 
one-instanton contribution, or first trans-series, yields an expansion of the form 



A 



J* 1 \z) = izV 2 e-^F 1>1 \^ + Y J F 1 , n+l z n l 2 ^ , (4.154) 

where z = g 2 . This is an important feature distinguishing matrix models and string theory from 
field theory and QM: the action of an instanton goes like 1/y/z, and not as \/z. Similarly, the 
perturbation series around the instanton sector is a series in powers of ^/z, and not a series in 
powers of z. We may now write 

oo 

J^°)(z) = ^2F g (t)z 9 . (4.155) 

Let us assume that there is a "resurgent" relation between the perturbative sector and the first 
trans-series. In analogy with the calculation in (2.133), we expect 



n_1 ' (4.156) 

*%i | A -W T{2g -p) + jr F w A- 2 ^+ n Y{2g - p - n) ) . 

71=1 / 



7T 



This leads to an asymptotic expansion of F g in powers of 1/g. Up to two loops we can write it 



9 

as 



F g — F(2g-P)F ltl 



F 12 A 
1 + ^ + 



(4.157) 



This gives a prediction for the large-order behavior of genus g free energies in the one-cut matrix 
model. 

In writing (4.157) we have implicitly assumed that there is a single instanton solution that 
contributes to the asymptotic behavior. In general there might be various instanton configura- 
tions in the system, with the same action in absolute value, and in this case F^ will denote the 
sum of all these contributions. It is also common to have complex instanton solutions which give 
complex conjugate contributions to F^\ as we saw in QM. In this case the asymptotic behavior 
of F g is again obtained by adding their contributions. If we write 

A=\A\e- i9 \ F 1>1 = \ t x\e w , (4.158) 

the leading asymptotics will read in this case 

F g ~ ] -> T(2g-f3)\fi\ co S ((2g- f3)9 A + e). (4.159) 

7T 
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Figure 28: The left-hand side shows the spectral curve in the one-cut phase of the cubic matrix model. 
The instanton action relevant in the double-scaling limit is obtained by calculating the B-period of the 
one-form y(x)dx, which goes from the filled cut A\ to the pinched point. The two-cut phase, in which 
the pinched point becomes a filled interval, is shown on the right hand side. The instanton action is still 
given by the B-pcriod integral. 



We have also assumed so far that the large order behavior is dominated by the instanton 
associated to eigenvalue tunneling. However, there is a subtlety concerning the large order 
behavior and related to the Gaussian part of the free energy. For small t, the behavior of the 
free energy is dominated by its Gaussian part, 

FJt) sa . B2g - „ , t->Q. (4.160) 
9K ' 2g{2g-2)t 2 ^ V ; 

Using the asymptotics (2.97), we can see that the Gaussian part indeed displays the large g 
behavior (4.157), with 

A = -27rt, (3 = 1. (4.161) 

However, this behavior is due to the universal Gaussian part of the free energy. For other 
values of t, the large order behavior is indeed controlled by the non-trivial instantons due to 
eigenvalue tunneling and discussed above. In the quartic matrix model, for example, there is a 
non-trivial instanton with action (4.99), and one sees from the explicit calculation in (4.84) that 
this instanton has 

P=\- (4-162) 

One can test the conjectural asymptotics (4.157) governed by this non-trivial instanton by con- 
sidering the normalized free energy obtained, i.e. by subtracting the Gaussian free energy. This 
was done numerically in [100, 10] by looking at the sequence of the first F g , and one can see that 
(4.157) is indeed satisfied. 

The appearance of two different instanton actions, corresponding to the "trivial" action 
(4.161) and the non-trivial instanton action associated to eigenvalue tunneling, has a natural 
geometric interpretation [52]. The instanton action describing eigenvalue tunneling in the one- 
cut phase can be written in terms of the spectral curve y(x) as in (4.86). This can be in turn 
written as a period integral of the natural meromorphic form y(x)dx along a B cycle which goes 
from the filled cut to the critical point: 

A B = ^^y(x)dx. (4.163) 

In Fig. 28 (left) we show the pinched curve describing the one-cut curve. The A% cycle corresponds 
to the filled cut, and the B cycle goes from A\ to the pinched cycle. The ^-period going around 



- 76 - 



the filled cut is just proportional to the 't Hooft parameter: 

A Al (t) = 2mt = \j A y(x)dx. (4.164) 

As we have just seen, both the A\ and the B periods give rise to instantons in the matrix model. 
The instanton corresponding to the A\ period is the "trivial" one and governs the asymptotics 
near t = 0, through the Gaussian part of the potential. In other regions of the t-plane, the 
large genus behavior will be controlled by 5-periods As{t) of the form (4.163). In general, the 
action controlling the large order behavior at a given point t will be proportional to the smallest 
period of the meromorphic form y(x)dx (in absolute value). The B-type periods As(t) vanish at 
critical values of the 't Hooft parameter, and the A\ period vanishes at t = 0, so in both cases 
the instanton action is given by a "vanishing cycle." 

This result can be generalized to two-cut phases, where the pinched point is now resolved 
into a second cut A2: the instanton action is still given by the S-cycle integral, now going from 
the first cycle A\ to the second cycle A2, see Fig. 28 (right). This instanton action controls the 
large order behavior of the free energies in the appropriate regions of moduli space. This was 
verified in [86] for the two-cut, cubic matrix model. 

This picture of instanton actions as periods of the spectral curve, pointed out in [52], seems 
to be the most general framework explaining the large order structure of a very general class of 
matrix models. For example, in the Chern-Simons matrix model for S 3 introduced in [91] and 
briefly reviewed in the Example 3.2, the instanton actions are given by [107] 

2vri (t + 2vrin) , n G Z. (4.165) 

For n = one recovers the action governing the Gaussian behavior. The instantons with n = ±1 
can be detected through the large order behavior of the genus g free energies, once the Gaussian 
part is subtracted [107]. Since the spectral curve describing this model admits a constant period, 
one can regard (4.165) as a linear combination of periods. A more detailed discussion of this 
general point of view on instanton actions can be found in [52] . 

4.7 Classical asymptotics and the Stokes phenomenon in matrix models 

In the previous subsections we have introduced the basic ingredients to study the large N asymp- 
totics of matrix models. Once a background is chosen, it leads to a formal perturbative series in 
g s , as in (4.15), and the rest of the sectors in (4.7) lead to a formal trans-series. 

We would now like to proceed to step two in the program sketched in the introduction, i.e. 
we want solve the problem of classical asymptotics in the case of matrix models at large N, or, 
equivalently, we want to find the formal expansions which provide the asymptotics in each region 
of the complex space of parameters, paying attention to the Stokes phenomenon. 

This turns out to be a non-trivial and rich problem. We will discuss some of its aspects by 
looking at a concrete example: we will consider the matrix model version of the Airy function. 
This model is defined by the partition function (4.4) where 7 = 71 is the contour shown in Fig. 2 
and the potential V(x) is given by (4.102). Without loss of generality, we will assume that g s is 
real and positive, and we will study how this partition function changes as we change t > (the 
't Hooft parameter) and the parameter k in the potential. 

For small t, as we discussed in (4.22), the saddle point structure should be the same as for 
the Airy function. Therefore, at small t we expect the following phase structure. For \k\ < 2n/3, 
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the dominant saddle is the one-cut configuration where all the eigenvalues sit near x L = — C 1 ^ 2 , 
i.e. 

(iV,0), (4.166) 

and 

Z ri (N,g s )~ Z(N,0). (4.167) 

The other saddles, where eigenvalues will tunnel to the critical point at x R , are not relevant since 
the integration path does not pass through x R . 

For 2ir/3 < \k\ < it the integration path gets deformed into the sum of the two steepest- 
descent paths, and the other instanton configurations will start contributing by eigenvalue tun- 
neling. However, at least at small t they should be exponentially suppressed. In other words, the 
dominant configuration to the matrix integral in this region should still be the boundary saddle 
(N,0), but there will be exponentially small corrections due to tunneling to the other critical 
point x R , and we will have 

Z 7l (N,g s ) = Z(N,0) i^ + fjO'Z^iN^j , (4.168) 

where Z&{N) is given by (4.118) and C is an appropriate constant. In contrast to the case 
of the Airy function, in the large N limit there is an infinite number of corrections involving 
the subleading saddle points. The above formula is the analogue of (2.73) for the "matrix Airy 
integral." Of course, from the point of view of classical asymptotics, this infinite series of expo- 
nentially small corrections is not taken into account, but it will be needed if we want to perform 
Borel resummations. Notice that the exponentially small corrections in (4.118) correspond to 
the eigenvalue tunneling depicted in Fig. 25, i.e. 

(N, 0) — >■ (N — £,£). (4.169) 

Now, a very interesting thing happens when we reach the anti-Stokes line k = tt. The 
real part of the action A(t) vanishes at t = 0, and it is actually negative for small t, therefore 
the subleading corrections are no longer negligeable. This is the familiar "return of the small 
exponential" along the anti-Stokes line. There is however a crucial difference with the case 
of the Airy function. The reason is that, in the large N limit, there is an infinite number of 
corrections that come to life on that line, namely all the sub- leading terms in (4.118). In fact, 
the same situation arises in non-linear ODEs, where this leads to the formation of singularities 
along anti-Stokes lines, see for example [40]. 

It is clear that in order to obtain a reasonable asymptotics we must somehow sum all the 
corrections. We proceed as follows. It is clear that the saddle (N, 0) is no longer a good starting 
point for the expansion. We have to find a new dominant saddle, which must necessarily have 
N% / 0. In order for this saddle to be stable, the effective potentials on both cuts have to be the 
same. This requires 

Re ( = 0, (4.170) 



9s 

where the derivative is w.r.t. the variable s defined in (4.119). This condition was first proposed 
in [80], and it is an equipotential condition which blocks the "flow" of eigenvalues. A way to 
derive this condition is to write the sum (4.7), in the case d = 2, as 

C^( ( ^) n e^F (n) + - (4171) 



n=0 
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where we denoted N\ = N — n, N 2 = n. Requiring this sum to have a saddle-point for n* leads 
to (4.170). In the large N limit, it has been proposed in [27] to replace the sum over n by a theta 
function. One then gets the following asymptotic expansion in g s , 

k mi>O gi>1 -^ 



(4.172) 

The derivatives of the free energies F g are again w.r.t. s. The theta function 6 M;i , with charac- 
teristics (/Lt, z/) is defined by 



Q/j,v(u,t) = ^ e {n+n-Ne)u e 7ri(n+ ti-Ne)r(n+fi-Ne) Q 2m(n+^)u ^ (4.173) 

and it is evaluated at 



9s 27T1 

In the above equation, we have denoted 



Fq. (4.174) 



e = f, C=| = e 2 -. (4.175) 

In standard matrix models we have /U = 0, although the term iVe might give an effective charac- 
teristic depending on the parity of N [27] . 

The asymptotics (4.172), discovered in [27], was interpreted in in [97] as a matrix model 
generalization of the oscillatory asymptotics along an anti-Stokes line. The singularities along 
anti-Stokes lines in nonlinear ODEs correspond here to zeros of the partition function, which 
are made possible due to the presence in (4.172) of the theta function at leading order: when 
the theta function vanishes, the partition function vanishes at leading order in g s . These zeros 
are then nothing but Lee- Yang zeros for the partition function of the matrix model, which are 
known to occur along anti-Stokes lines [78, 107]. 

Notice that the term involving the theta function is not analytic in N, in the same way that 
the asymptotics (2.76) is not analytic at x = oo (even when expressed in terms of z = X s / 2 ). We 
have, for the free energy, the expansion 

F = g; 2 F (N*) + F 1 (N*) + log 9^ + • • • , (4.176) 

so the oscillatory behavior is already present at next-to-leading order. Notice that this asymp- 
totics has contributions which are not present in the standard large N expansion (4.15). Let us 
now come back to the cubic matrix model along k = ir. We have to find a new saddle satisfying 
(4.170). For arg(ft) = tt, C 1 / 2 is purely imaginary, and it is easy to see from the structure of the 
genus zero free energy that its real part is symmetric in t±, t 2 . Therefore, 

h = h (4.177) 

solves the saddle-point equation 

f)F 

Re— - = 0, (4.178) 

OS 
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at least for small t. We conclude that, on the anti-Stokes line and for t small enough, the 
partition function of the cubic matrix model is given by an expansion of the form (4.172). Since 
C\ = C2 = 1 (this follows from the fact that the path 71 is deformed into a sum of two steepest- 
descent paths), the characteristics of the theta function are \x = v = 0. 




Figure 29: The phase diagram of the cubic matrix model with potential (4.102), as a function of 2ir/3 < 
K < n and real t = g s N. The parameter t is also supposed to be small. The Stokes and anti-Stokes lines 
of the N = 1 case, which occur for n = 2tt/3 and k = 7r, respectively, determine to a large extent the 
phase diagram for small t. On the anti-Stokes line, the saddle value of t 2 is t/2, and the classical l/N 
asymptotics requires corrections beyond the genus expansion. 

The above analysis is only valid, strictly speaking, at t ~ 0. When t is small but nonzero, 
the real part of the instanton action vanishes at a finite t c {n) for each value of 2tt/3 < k < it. 
Therefore, the transition to a new saddle will occur along a line in the t — k plane defined by 

Re(A(t c (K),K)) = 0. (4.179) 

We display the phase diagram of the model in Fig. 29. For a generic point along the curve t c (n) 
separating the two phases, there is a phase transition in which N% increases smoothly away from 
zero (this has been verified very explicitly in another, related model in [97]). This transition is 
then of the type "birth-of-a-cut" studied in for example [55]. Notice that the smooth transition 
becomes discontinuous as t — > i.e. it becomes the jump in asymptotics along the anti-Stokes 
line. For t > 0, the discontinuity is smoothed out. We then have a deformed Stokes phenomenon 
as we turn on the 't Hooft parameter. Notice that the phase diagram that we have just sketched 
is valid for small t. For arbitrary t there is a very rich phase structure first described in [44] and 
which can be formalized in terms of so-called Boutroux curves, see [23]. A recent analysis of the 
phase structure of the cubic matrix model can be found in [7]. 

Based on this detailed example, we can extract some more general conclusions for the study 
of asymptotics and non-perturbative effects in matrix models: 

1. As we change the parameters of the model, the partition function (4.7) is dominated by 
different large N saddle points in different "phases". There are phase transitions among 
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these different phases. The transitions occur as in asymptotic analysis: as we move in 
parameter space, large N saddles which were exponentially suppressed are no longer so and 
start contributing to the asymptotics. These transitions, triggered by large N instantons, 
are called large N phase transitions, and occur along generalized anti-Stokes lines. 

2. The Stokes phenomenon is "smoothed out" at finite 't Hooft parameter, even in the sense 
of classical asymptotics. 

3. The 1/N expansion breaks down along anti-Stokes lines. In the Hermitian case, the classical 
asymptotics involves theta functions with a non-analytic dependence on N. 

A more detailed study of these phenomena, together with an analysis of the so-called lens space 
matrix model, can be found in [97]. 

5. Two applications in String Theory 

It is impossible to provide here a detailed overview of non-perturbative effects in string theory. 
We will content ourselves with a discussion of some string theory models where one can use large 
N matrix models to describe non-perturbative effects. 

From the point of view discussed in these lectures, it is useful to focus on string free energies 
and to regard string theory as a quantum mechanical system with two different Planck constants. 
The worldsheet Planck constant is given by the square of the string length, 

Ks = £ 2 S . (5.1) 

The free energy at genus g, F g (t), depends on a set of target moduli t (this could be a compact- 
ification radius or a Kahler parameter), and the perturbative regime of the worldsheet theory 
corresponds to the regime in which these moduli are very large as compared to the string length. 
Non-perturbative effects with respect to this Planck constant are of the form 

exp (-A ws (t)/£ 2 S ) , (5.2) 

where A ws (t) is a worldsheet instanton action. These effects are usually realized in string theory 
as instantons of the non-linear sigma model underlying perturbative string theory. 

There is however a second, spacetime Planck constant called the string coupling constant, 

tht = 9s- (5.3) 

The contribution of a genus g worldsheets to the free energy is weigthed by the factor g s , and 
the perturbative free energy is given by a genus expansion of the form 

oo 

F(°Ht,g s ) = J2F 9 (t)9 2 s 9 - 2 . (5.4) 

This has the same structure that the large N expansion of a gauge theory or a matrix model, 
compare to (4.15), and it is believed that, generically, the genus g free energies grow as (2g)\, 
just as in large iV theories [71, 112]. We should then expect non-perturbative effects w.r.t. this 
second Planck constant, of the form 

exp{-A st (t)/g s ) (5.5) 
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where A st (t) is the action of a spacetime instanton. These effects were first found experimentally 
in non-critical string theory [44, 112], and later on it was proposed by Polchinski that they are 
realized by D-branes [108]. 

In general, the calculation of these non-perturbative effects directly in string theory is not 
easy. However, when the string theory model has a large N gauge theory dual, one can in 
principle compute them in the gauge theory. Worldsheet instantons appear as exponentially 
small corrections in the 't Hooft expansion at strong 't Hooft coupling, while spacetime instantons 
should correspond to large TV instantons. In some cases, the large ./V dual further reduces to a 
matrix model, and the techniques described in these lectures can then be applied to calculate 
non-perturbative effects in string theory. In particular, the free energy of non-critical strings is 
described by non- linear ODEs, and non-perturbative effects in these theories can be analyzed 
with the trans-series formalism presented in section 2 (see also [35] for recent developments 
on non-critical strings). We will now describe two examples where the ideas developed in these 
lectures lead to non-perturbative results in string theory, and which go beyond non-critical string 
theory. 

5.1 A toy model: Trans-series and Hurwitz theory 

Hurwitz theory can be regarded as a toy model for string theory (or, more precisely, for topological 
string theory). In Hurwitz theory one considers maps from a Riemann surface or string, to 
another Riemann surface, and "counts" these maps in an appropriate form, by using Hurwitz 
numbers. Here, for simplicity, we will assume that the target is P 1 , i.e. a two-sphere. The only 
parameters in Hurwitz theory are the string coupling constant gjj and a target parameter tjj 
which can be regarded as measuring the "size" of the target. The partition function of this model 
is of the form 

where Q = e~ tH . The dependence on gjj is the standard one for a sum over genera in string 
theory. The quantity d {l d ) is a simple Hurwitz number counting degree d covering maps of 
P 1 , with simple branch points only, and by Riemann surfaces of genus g (see for example [100] 
for an explicit expression for this number in terms of representation theory data). Here, in the 
partition function, we consider maps from generally disconnected worldsheets. The free energy 
F = log Z describes connected, simple Hurwitz numbers H^ d {l d )' , 

F(g H ,t H ) = ^9 2 H g ' 2 F g (Q), (5.7) 
<?>o 

where 

- H pl (l d )* 

This theory is in fact a string theory in disguise. It can be realized as a special limit of topological 
string theory on certain toric Calabi-Yau manifolds, see for example [30, 33] for detailed deriva- 
tions. It was also conjectured in [91] and proved in [56] that Hurwitz theory can be described 
in terms of matrix integrals, and this in turn was used in [100] to compute instanton effects, by 
using the techniques reviewed in subsection 4.5.1. As shown in [94] there is however another 
way to understand non-perturbative effects in this theory, by using difference equations as in 
subsection 4.5.3. 
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As shown in [105], the free energy of Hurwitz theory satisfies a difference equation of the 
Toda type, 

exp(F(t H + g H ) + F(t H - g H ) - 2F(t H )) = gW H d? H F(t H , 9h)- (5.9) 

As we did in subsection 4.5.3, we can try to solve this equation with a trans-series ansatz for the 
free energy of the form (4.144). Doing this one immediately obtains the following equation for 
the one-instanton amplitude, 

exp(A te F(°)(^)) ^ 9H F^(t H ) = gWdlF^, (5.10) 

where we have written 

A h /(t) = /(* + h) + f(t -h)- 2f(t) (5.11) 

to denote the discrete Laplace operator with step h. The first term in the expansion of (5.10) in 
powers of gu gives an equation for A'^h), 

cosh(A'(* H )) - l] = exp{t H - dlF^(t H )} (A'(t H )) 2 . (5.12) 

It is straightforward to use (5.9) to derive a one-parameter, trans-series solution for the free energy 
of the Hurwitz model. In addition, one can check that the function A(tjj) defined implicitly by 
(5.12) coincides with the instanton action computed in [100] from the matrix model realization 
of [91, 56]. This instanton action controls the large genus behavior of the free energies F g {Q) of 
Hurwitz theory. We conclude that, in this toy model of string theory, the perturbative sector 
obtained by considering Hurwitz coverings can be generalized to a general trans-series. However, 
the geometric or physical meaning of this trans-series is not yet known. 

The Hurwitz model that we have discussed can be generalized to a family of topological 
string theories on non-compact Calabi-Yau manifolds (the so-called "local curves") where non- 
perturbative effects can be also computed by using matrix model techniques, see [93, 100, 95]. 

5.2 Stringy instantons in ABJM theory 

We will now describe an application in physical superstring theory, namely ABJM theory [3] and 
its type IIA dual. 

ABJM theory [3, 12] is a quiver Chern-Simons-matter theory in three dimensions with gauge 
group U{N)k x U(N)-k an d N = 6 supersymmetry. The Chern-Simons actions for the gauge 
groups are of the form (3.51) and they have couplings k and —k, respectively. The theory contains 
as well four bosonic fields Cj, / = 1, ■ ■ ■ , 4, in the bifundamental representation of the gauge 
group. The 't Hooft parameter of this theory is defined as 

TV 

A = -. (5.13) 

In [82] it was shown, through a beautiful application of localization techniques, that the partition 
function of ABJM theory on the three-sphere can be computed by a matrix model (see [96] for 
a pedagogical review). This matrix model is given by 



I /•^dwd^n i <,smh 2 (^/ 1 )sinh 2 (^) 

w«, 9 ,> -^/n^ n ,^m ^ 
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where the coupling g s is related to the Chern-Simons coupling k of ABJM theory as 

9s = (5.15) 

The reduction of the partititon function of the theory to a matrix integral makes it possible to 
address some of the questions considered in these lectures in a very concrete way. For example, 
one can wonder, as we did in subsection 3.2 for Chern-Simons theory, what is the large order 
behavior of perturbation theory for the partition function or the free energy. For fixed N = 2, 
it has been shown in [109] that the perturbative expansion in powers of g s grows factorially but 
it is Borel summable. A different, more complicated question concerns the behavior of the free 
energy in the large A expansion, i.e.: what is the behavior of the genus g free energies? What 
are the non-perturbative effects at large .A? These issues were addressed in [51, 52], which we 
now review. 

The above matrix integral turns out to be a close cousin of the matrix integral for CS theory 
mentioned in (3.75). Although in principle it is not of the form (4.4), it can be put in such a 
form by an appropriate change of variables. Moreover, all the notions introduced in the previous 
section to analyze the large N limit of conventional matrix models -such as resolvent, spectral 
curve, and the like- can be generalized to the ABJM matrix model. It turns out that, at large 
A, the ABJM matrix model is a two-cut matrix model. In terms of the variables 

Y = e v , X = e x , (5.16) 

the spectral curve is given by the equation [2, 74, 99, 51] 

X 2 

Y + — - X 2 + \kX -1 = Q. (5.17) 

The Riemann surface of (5.17) can be represented by two A-planes glued along the cuts [1/a, a] 
and [—6, —1/6]. The position of the endpoints can be determined from 

a + - + b+\=A, a+ l -b-\ = 2\K. (5.18) 
a b a b 

The variable k can be regarded as the complex modulus of the spectral curve. As in (4.49), the 
't Hooft parameter can be obtained as a period of the spectral curve and it is related to k by [99] 

k „ (\ 1 1 3 k 2 



AW = S^2<2-2 ;1 -2 ; -I<!j- (5 ' 19) 
The free energy of the ABJM matrix model has a 1/A expansion of the form 

oo 

F{\g s ) = Y j g 2 r 2 F g {\). (5.20) 

9=0 

The genus zero free energy can be obtained from (4.50) as [51, 96] 

^°-I G 3,3^ 0; 0j _i -Tq) + ^T 3jP2 \2> 2' 2' X ' 2' "16 J' (5 - 21) 

2 3 

wheere G 3 ' 3 is a Meijer function. Near A = this planar free energy has a logarithmic singularity, 
but this is due to the Gaussian part of the free energy, and once this part is subtracted (as in 
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(4.20), with d = 2 since there are two cuts) we obtain an analytic function at A = 0, as expected. 
There are singularities in the complex plane of the 't Hooft parameter, signaling the finite radius 
of convergence of the series. The location of the singularities can be found explicitly from the 
above expressions. Indeed, the singularities correspond to the branch point of the hypergeometric 
functions involved, at 

k = ±4i, (5.22) 

which correspond to 

2iif 

A = T — o - , (5-23) 

where K is Catalan's constant. For reasons explained in [51], these points are called the conifold 
points. 

In [51] it was also shown that, for g > 1, the free energies can be written in terms of 
quasi-modular forms of the modular parameter of the elliptic curve (5.17). 

K^fl ( 4) 

For g = I, one simply has 

Fi = -Iog^(r), (5.25) 

where i] is the usual Dedekind eta function. For g > 2, the F g can be written in terms of the 
quasi- modular forms E2 (the standard Eisenstein series), b and d, where 

b = 4(r), d = #j(r), (5.26) 

are standard Jacobi theta functions. More precisely, we have the general structure 

1 3g-3 

F ^ = Tr^i E El(r)pi 9 \b,d), g > 2, (5.27) 

where p^\b, d) are polynomials in 6, d of modular weight 6g — 6 — 2k. In [51] a recursive procedure 
was described which gives all the genus g free energies unambiguously. The genus g free energies 
Fg(X) obtained in this way are exact interpolating functions of the 't Hooft parameter, and they 
can be studied in various regimes. When A — > they reproduce the perturbation theory of the 
matrix model (5.14) around the Gaussian point A = 0, where they behave, as expected, as two 
copies of the Gaussian matrix model, 

FaW = - g(2ff-2) (27riA)2 ~ 2g + ° {X) - (5 - 28) 

The genus g free energies can be also studied in the strong coupling regime A — > 00, or, equiva- 
lently, at k — > 00. In this regime it is more convenient to use the shifted variable 

24 ~ 2w 2 

One finds the following structure. For Fq and F\ one has, at strong coupling, 



i = i -^7 = 1 ^+0(^ 2 ), (5.29) 



Fl = \ l0gK ~\ l0g 



-27TV2A 



21ogK 


+ ( 
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(5.30) 
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For g > 2 one has 



(5.32) 



where 

A9- l \B 2g B 2g _ 2 \ 
g(2g-2)(2g-2)\ 

can be interpreted as the contribution from constant maps to the free energy, and 

f 9 (x) = Y,4 9)x29 ~ 3+j ( 5 - 33 ) 

3=0 

is a polynomial 1 . The leading, strong coupling behavior is then given by 

Fg(X) -c g ~ Ai~ 9 , A -»• oo, g > 2. (5.34) 



We can now ask what is the large order behavior of F g (X) and what are the possible instanton 
configurations underlying this behavior. First of all, we notice that the A-independent part c g of 
F g has the large genus behavior, 

c 9 ~ -(27r 2 r^T(2g-l), (5.35) 

and corresponds to a constant instanton action governing this constant map contribution [107], 

A m = 2tt 2 . (5.36) 

In order to identify other instanton actions contributing to the asymptotics, one can look at 
vanishing periods of the meromorphic form y(x)dx, as suggested by the discussion in subsection 
4.6. There are two obvious vanishing periods, namely the "Gaussian" period vanishing at A = 0, 

A w {k) = -4tt 2 A(k). (5.37) 

There is another period which vanishes at the conifold point k = ±4i, given by 



A (k) - + 4^ 2 A ± vr 2 - — G 2 ' 3 I 2' 2' 5 



^ 2 * _2 



iri) ±vr 2 . (5.38) 



Finally, there is a linear combination of both instanton actions, 

A s {k) = A w {k) + A c (k). (5.39) 

It was shown in [52] that these three instanton actions control the large order behavior of F g (X) — 
c g in different regions of the complex plane of the 't Hooft parameter. Near the Gaussian or 
weakly coupled point, it is (5.37) which controls the large order behavior. Near the conifold 
point (5.23), it is (5.38) which dominates. Finally, in the strong coupling region A » 1, the 
large order behavior is controlled by (5.39). Notice that we can regard the behavior of (5.35) as 
controlled by a constant period, which always exists in these models. 



lr The interpolating functions F g computed in [51] did not include the constant map contribution. This was 
corrected in [76]. 
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We can now discuss Borel summability. In the physical ABJM theory, A is real and g s is 
purely imaginary. The expansion (5.20) should be written in terms of the real coupling constant 
2ir/k, i.e. as 

F ( A > fc ) = E t (-ir 1 ^)- (5.40) 

We get an extra (— l) 9 " 1 sign at each genus. Equivalently, this leads to an extra — i factor in 
the instanton actions computed above. We can now ask whether this factorially divergent series 
is Borel summable or not. At strong coupling A > 1, the behavior of the genus g free energy 
(— l) 9 ~ 1 F g (X) is dominated by the constant period —\A m = — 2i7r 2 , which is purely imaginary. 
We then obtain a Borel summable series. Even after subtracting the constant map contribution, 
we obtain a Borel-summable series at strong coupling. Indeed, the strong coupling action (5.39), 
which controls the asymptotics in this regime, is complex: 

Im(-L4 S (A)) = vr 2 , (5.41) 

and for large A we have, 

-L4 S (A) = 2iT 2 V2X + TT 2 i + (e~ 2w ^ , A > 1. (5.42) 

Interestingly, it was shown in [52] that the first, leading term of this action at strong coupling 
coincides with the action of a D2-brane wrapping MP 3 inside CP 3 . This seems to indicate that 
this instanton of the matrix model corresponds indeed to a D2-brane in the dual type IIA string 
theory, as expected from general arguments in superstring theory [108]. 



6. Concluding remarks 

In these lectures I have tried to provide a pedagogical introduction to some aspects of non- 
perturbative effects in quantum theory. I have focused on topics which are not covered in detail 
in the classic reviews on the subject, like large N gauge theories and matrix models, with a view 
towards applications to string theory. It is clear that, in many respects, this topic is still in its 
infancy. Even in the context of large N matrix models, where many explicit results are available, 
there are still some important open problems. For example, there are no detailed results on the 
large order asymptotics for the genus g free energies in multi-cut phases, i.e. there is no analogue 
of equation (4.156) for these phases; only the instanton action A has been identified in some 
cases, as discussed in subsection 4.6. The ideas of resurgence, which have been very powerful in 
the context of ODEs and Quantum Mechanics, have scarcely been used in QFT (see however [11] 
for recent work in this direction). In the case of string theory, these ideas have been explored in 
[94, 10], but mostly in the case of non-critical strings and toy models with large N duals, and 
much remains to be done. I only hope that these lectures will be useful for future research along 
these directions. 
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